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1.0 SUMMARY 


An analytical model for aircraft interior noise prediction is 
considered in this report. The model can be used to predict the 
sound levels inside an airplane cabin caused by the rotation of 
a propeller (of any design) alongside. The fuselage model is 
that of a cylinder with a structurally integral floor. The 
cabin sidewall and the floor are stiffened by ring frames and 
stringers or floor beams of arbitrary configurations. The cabin 
interior is covered with a trim (i.e., layers of insulation and 
septa with a lining) to increase the sidewall sound isolation 
and provide absorption in the cabin. 

The results are the culmination of a three phase program 
sponsored by NASA Langley Research Center. In Phase I the basic 
analytical modeling of the transmission problem (interaction of 
the structure with the exterior and interior acoustic fields) 
was undertaken and preliminary validation studies were completed 
using an unpressurized, unstiffened cylinder as a test article. 
Results of that work are presented in Reference [3]. In Phase 
II, the general aircraft interior noise model was developed and 
preliminary work on the laying out of the basic master computer 
program began. Validation studies were conducted using more 
advanced test articles (one being a stiffened cylinder with a 
floor partition and interior trim). Results of that work are 
found in References [4] and [7]. In Phase III, the analytical 
models and the software were completed (including the propeller 
excitation work). Validation studies using a scale model fuse- 
lage excited by a propeller were undertaken and the documenta- 
tion of the finalized model and software package was completed. 

The present model is believed to be the only one in existence 
that can be used to calculate the interior sound levels using as 
input data, the precise propeller noise signature over the 
fuselage. 
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2 . 0 INTRODUCTION 


This report presents the details of a basic airplane interior 
noise model. The elements of this model include a fuselage and 
a propeller (Figures 1 and 2) . The fuselage consists of a cy- 
linder stiffened by ring frames and stringers, and a floor that 
is structurally an integral part of the fuselage. The cabin 
space is the volume above the floor. The interior surface of 
the cabin (sidewall) is finished out with a trim consisting of 
insulation covered with a lining. The propeller rotates about an 
axis parallel to the centerline of the fuselage. The model can 
be used to predict the sound levels in the cabin space for each 
of the various harmonics of the propeller. 

The excitation of the exterior of the fuselage is obtained using 
a propeller noise prediction model developed by NASA Langley. 

The present model works with the pressure time histories (signa- 
tures) as defined over the fuselage at a number of closely 
spaced points on a grid that lies in the fuselage skin. The 
pressure signatures are Fourier analyzed to define the ampli- 
tudes and phases of each of the harmonics of the propeller tones 
(at each location on the grid). The cross power spectral den- 
sity function for each harmonic, for all grid point pairs (a 

\ 

delta function in the frequency domain) is used to compute the 
values of the generalized forces for each structural mode of the 
fuselage. 

The fuselage structural modes are developed for the case of a 
stiffened cylinder with a floor partition. The structural modes 
are described by their eigenvalues (resonance frequencies), 
eigenvectors (mode shapes), and loss factors. The mode shapes 
include not only the cylinder wall normal displacement (w com- 
ponent) but also the normal displacement of the floor, and the 
in-plane axial and circumferential displacements (u and v 
components) of cylinder and floor as well. The loss factors of 
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FIGURE 2. PROPELLER AND FUSELAGE SURFACE POINT GEOMETRY 


the bare fuselage are input and must come from measurements. 

When trim is installed on the sidewall, the structural losses 
increase due to the trim's presence against the sidewall and 
this is computed for the particular trim installation. 

The displacement of the trim lining induced by the fuselage skin 
vibration is determined using a transfer matrix which relates 
the pressures on the inside of the skin and the inside of the 
lining to the displacement of the skin and of the lining. The 
transfer matrix contains all of the physical properties of the 
insulation and lining required for the calculation. The wave 
impedance of the insulation and complex acoustic wavenumber are 
input as physical parameters to describe the insulation. 

The coupling of the lining to the interior acoustic field is 
calculated for each acoustic and structural mode. The acoustic 
modes are defined by their resonance frequencies, mode shapes, 
and loss factors. The acoustic loss factors must be input for a 
bare fuselage but are calculated when a cabin trim is installed, 
from elements of the trim transfer matrix. 

The model allows for the calculation of the space average mean 
square pressure in the cabin for each propeller harmonic (up to 
a maximum of ten (10) harmonics). 

2.1 Report Organization 

This report considers the analytical derivations, experimental 
tests, and validation studies. The analytical derivations are 
presented in Section 3 and Appendices A through D. Appendix E 
is devoted to test comparisons and to the determination of the 
quality of the predictions. Appendix F is a list of symbols 
used. 

Early in Section 3, general solutions are given for the basic 
sound transmission problems of concern. Problems of tone 
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transmission (propeller) and noise transmission (reverberant 
exterior field) are considered. Solutions that are to be used 
in the low and high frequency regimes are presented (the low 
frequency formulations apply until the acoustic modal density of 
the cabin space is equivalent to 10 to 15 modes in every one- 
third octave band). Beginning with Section 3.4, the various 
terms appearing in the general solutions are evaluated, literal- 
ly specialized to create the desired interior noise model. 

These terms include the generalized forces for propeller noise 
excitation, interior coupling factors, joint acceptances, 
resonance frequencies and loss factors. Information needed to 
complete these calculations are derived in the first four 
appendices . 

The model validation undertaken in Phase III is considered in 
Appendix E. The experiment and the data acquired are discussed 
and a statistical comparison of predictions and measurements is 
presented. 


2.2 Program Management 


The work was accomplished in joint effort by BBN/Los Angeles and 
NASA Langley Research Center. The experimental work was done at 
NASA Langley by C. M. Willis and W. H. Mayes. Mr. Mayes acted 
as LaRC technical representative of the contracting officer 
(TRCO). L. D. Pope served as BBN program manager. 
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3.0 ELEMENTS OF THE PROPELLER AIRCRAFT INTERIOR 
NOISE MODEL 

A detailed description of the propeller aircraft interior 
noise prediction model is given in the following paragraphs. 
It is the intention of the authors to present all of the 
background needed to understand the equations and their 
origins. This is done partly through the use of a series of 
Appendices that summarize the primary results from two 
previous phases and the now concluding phase of work that 
have led to the present model. References to some published 
papers and books are also required. 

A concern in this report is to define a reasonable level at 
which to begin the technical presentation. This model is 
based on the general solutions of sound transmission prob- 
lems for cases of tonal and broadband noise excitations. 
Parts of the report deal with the theory of fluid-structure 
interaction problems with direction at the development of a 
sufficiently general solution as needed for the present 
purpose. Because the propeller is an integral part of the 
model the field of "aeroacoustics" is inherently involved. 
Due to the complexity of the fuselage models and the geomet- 
ry of the cabin space, the structural dynamics of stiffened 
shell structures must be considered and numerical procedures 
for calculation of the acoustic characteristics of complex 
spaces devised. The presence of sidewall trim (insulation 
and lining) leads to the use of the properties of porous 
materials . 

The ultimate user of this model will not necessarily be 
knowledgeable in all of these areas. For instance, one well 
versed in the fluid-structure intraction problem may have 
little or no experience in propeller noise. Keeping this in 
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mind, an effort has been made to ease the transition between 
the various disciplines and the introduction to each topic 
is kept at what might be considered an intermediate level. 

3.1 General Solution for the Sound Transmission Problem 


The fundamental goal of this model is the prediction of the 
sound level inside an airplane cabin due to noise on the 
exterior caused by the rotation of a propeller. The 
propeller noise of concern is at discrete frequencies and is 
not really a noise at all, rather a series of tones. The 
lowest frequency tone is at the blade passage frequency 
(BPP) and the other tones are its harmonics, that is, they 
occur at frequencies that are integral multiples of the BPF. 
As a secondary goal of this model, predictions for the cases 
of an arbitrary exterior harmonic (tonal) field and for a 
reverberant exterior (noise) field are sought. Sound 
transmission through the fuselage sidewall and the trim 
insulation and lining is of concern. 


The basic method used herein to solve the sound transmission 
problem is that of a power balance. 


W in W diss 


( 1 ) 


The band-limited net time-averaged power, VT- n , flowing into 
the cabin must equal the net time-averaged power, 
dissipated on the cabin walls. By expressing the inflowing 
power in terms of the exterior exciting pressure and the 
dissipated power in terms of the interior response pressure 
and equating the expressions according to Eq.(l), one can 
solve for the interior pressure. 
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A primary feature of the power flow approach adopted is that 
integration of spectral components has been performed 
analytically to achieve band-limited levels, i.e.. 


w. = 

m 

L 

Re | 

[ 0 ->] 

dm 

diss 

4 

Am 

Re 

D&>] 

| dm 


(2) 


W ra d(“) is the complex spectral density of the power radiated by 
the structure into the interior acoustic space and W^g(w) is 
the spectral density of power absorbed on the inner wall ovf the 
space from the interior acoustic field. In the case of tone 
transmission, these quantities contain delta functions in the 
frequency domain. Eqs. (2) have applicability at all frequen- 
cies and can be applied in discrete modal representations even 
when there are no acoustic or structural modes resonant in the 
band Aoj. 


Power flow into each individual acoustic mode is computed, 
whether the mode is resonant in the band or not. If the mode is 
resonant in the band, one can speak of the "resonant response" 
of the mode; if the mode is resonant outside the band, the 
"nonresonant response" of the mode in the hand Ac. Similarly 
power flowing out of the acoustic field to the inside wall is 
computed mode-by-mode. By equating the power into a mode to the 
power out, the space-average mean-square modal pressure for mode 
n, limited to band Aco, is obtained, i.e., < p^> , , where s and t 

indicate the space and time averages, respectively. The band- 
limited, space-average mean-square pressure in the interior is 
the sum 


<pt>„ * = 5>P‘> 


i s , t 


n 


n s ,t 


( 3 ) 


9 



Here n is the complete set of acoustic modes, i.e., those 
resonant inside Aw (denoted by the symbolism, nsAm ), as well as 
those resonant above (n>Aw ) and below Aw (n<Aw ). It is noted 
here that the letter "n" always denotes acoustic modes and the 
letter "r" structural modes, usually being used as subscripts, 
superscripts, or indexes. 

Power Plow 


The concern is with a closed volume V, into which sound energy 
is propagating. The excitation of the space occurs at its 
boundary with the enclosing structure which vibrates due to an 
exterior source. For the present, the excitation can be consi- 
dered to be a broadband random acoustic field. For the specific 
case where the noise reduction is desired, the exterior field is 
taken to be diffuse (reverberant). 

The fundamental equation for the inflowing power comes from 
results developed in Refs. [1] and [2]. Specifically, use is 
made of Eq.(3) of Ref. [2]: 


W 


pc 


2 A 4 


in 


V 


S P bl U) 


j;<«) 

M 2 w 4 
r r 


t nV 2(n - 1 ' ) 


w„ 


x f 

1 { [1- 
Aw 


w 


n 
3 * 


n 


co 


(co 2 /w 2 )] 2 + n 2 H[l-(co 2 / C 0 2 J ] 2 +n£} 


(4) 


is the band-limited time averaged power which is being 

accepted by the acoustic field inside the enclosure, that is, 

the power received over frequencies lying in band Aw . The band 

has a width Aw = c w where co is the center frequency of the 

w 

band. c^ is a constant percentage which defines the bandwidth 
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(c w = 0.232 for one-third octaves). The quantity S Dbl (w) 

= <p^ 1 >/c w co is the power spectral density of the average 
exciting "blocked pressure" which acts on the transmitting 
structure. < Pbl > is the average mean-square pressure as would 
be measured in band Am obtained with a microphone sampling the 
surface pressure on the (hypothetical) immobilized (blocked) 
structure. The transmitting structure has an area A exposed to 
volume V. The average surface mass of the structure is m. The 
structure has a mass law sound transmission coefficient = 

(2pc 0 /wm) 2 , where p and c Q are the density and sound speed of 
air. 

The normalization of structural modes is taken such that the 
modal mass is 

M r = /x > 

V — — 

where \p (x) is the structural mode shape (mode r) at location x. 
x ranges over the entire structure, both transmitting and 
non-transmitting surfaces. Thus, if the rth mode of the 
transmitting structure (fuselage) is restrained in any way, say 
perhaps by the cabin floor, M r will take on a large value that 
will suppress transmission by that particular mode. The 
normalization of the acoustic modes is such that 



where <)> (£) is the value of the cabin acoustic mode shape (mode 
n) at location f . 

The coupling of the exterior field to the rth mode of the 
transmitting structure is given in terms of the joint acceptance 
j r (w). In the case where the exterior field is a noise it is 
assumed that j£( w ) changes little across Aw, i.e. by evaluating 
j£(w) at the center frequency of the band, a good approximation 
is obtained all across Aw. For this reason Aw must not be too 
wide, typically one-third octave or less. If the exterior field 
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is reverberant, j^(m) takes the appropriate form for that field. 
The same would be true for a progressive wave field. 


The coupling of a structural mode to an acoustic mode of the 
volume is given by the term 

f ' (n,r) = ~ 

In the above, x represents a point on the surface of the 
transmitting structure. Note that f'(n,r) is non-dimensional 
(hence the prime as a reminder) and that f ' (n,r) is always 
squared. f' 2 (n,r) is a positive number less than unity. The 
power flow depends on the location of the resonance frequencies 
m n and w r , i.e. relative to each other and to the band Am, and 
upon the acoustic and structural loss factors b n and b r . 


k 


4) n (“)^ r (x)dx 


Consider now the right-hand side of Eq . (1). W diss is the 
band-limited, time-averaged power that flows out of the interior 
acoustic field to the inner wall of the enclosure. This power 
flow is given by Eqs. (6) and (8) of Section III of [2] which 
reduce to 


W 


diss 


pc; 


n n M n 

0 ) 


<p 2 > 

n s,t 


n 


(5) 


where <p > is the interior space-average mean-square 
n s ^ u 

in band Aw attributable to the nth acoustic mode. The 
ment of Eq . (5) is presented in the appendix of [2]. 


pressure 

develop- 


In the case of a random exterior field, in which S_., ( co ) varies 

Pbl 

slowly in frequency, j ^, ( ) can be brought outside the integral 
as done in Eq . (4). However, for the general circumstance where 
tones are present, j£(cu) must be included under the integral 
sign, and replaced by its defining relation 


i 2 r M = 


//C pb i(x,x ;ii))4> r (x)ip r (x )dxdx 

^ S pbl (w) 
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where C (xlx'toj) is the cospectral density function of the 
pbl 

blocked exterior pressure field. 

Also, in order to get Eq . (A) in the desired form for the 
present model it must be modified in two other respects. First, 
to include fuselage sidewall trim, Eq.(4) is modified to include 
factor t , i.e. the trim transmission coefficient, to account 
for transmission through the trim, and the structural loss 
factor ri is replaced by to account for the added 
damping of the sidewall when trim is present. Second, proper 
account of the influence of the internal radiation damping of 
structural modes leads to the introduction of another term, n". 
The analytical developments for t, and n' are given in Appendix 

n M "k 

A; ,: r will be discussed in a subsequent section. After some 
rearrangement and the inclusion of the above results, Eq . (4) 

can be written in the following general form 



A bar over f'(n,r) indicates trim factor has been introduced. 
Eq . (6) is basically the same as Eq . (21) of Ref. [3], except 
modified to include trim and the proper influence of the 
internal radiation damping of resonant structural modes closely 

coupled in frequency to resonant acoustic modes. 
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The solution of a sound transmission problem is obtained by 

setting the nth term of Eq . (5) equal to the nth term of Eq.(6), 

solving for the individual modal pressures <p*> . one-by-one, 

n s y X/ 

and by adding the results according to Eq. (3), to obtain the 
interior level. 

3.2 Transmission of a Tone 


The fundamental model of concern here is one having to do with 
tone transmission. The propeller noise is, as stated, at 
discrete frequencies, each tone being a harmonic of the BPP. 
Let the various harmonics be , . . . , o> jj where 

“l = 2 ir x BPF 
oj 2 = 2 w i 
a) H = H^i , 

where H is the harmonic index. For a tone at frequency 
the co-spectral density function is defined by 

Cpbi ( x | x * ) 6 (o>-g> h ) = Cp bl (x|x';co) , 

C Pbl^ X l X 5 — ^ Cp^^(x|x T jco) j 0 < w < 00 , 



where Rpbi(5c |x' > T ) is the average cross correlation of the 
pressure over the blocked (immobile) fuselage, given by 

T 

R Pbi (x|x';x) = ^ J Re |~Pbi (x ,t )J Re[j> bl (x ' , t+x )J dt. 

"2 
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Inserting the first of the relations into Eq . (6) leads to the 
expression for the inflowing power to the cabin for the 
propeller tone at frequency ojh- 



H H — 

In this result, T ML and T t (in f'(n,r)) are evaluated at o» H . 

The space average mean square modal pressure is computed on a 
mode-by-mode basis using Eqs. (1), (5), and (7) (where w = <*> H ) : 



The interior space average level is then obtained from Eq . (3). 
Note that since 
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where <P^ is the mean square modal amplitude, the modal ampli- 
tude is determinable from 


<p 2 >^ = e <p 2 > H , 
n n ^n s,t 


3 


and theoretically, the mean square pressure at every interior 
point can be obtained with 


H 


<p 2 (C ,o))> t 


■ Dn<P’>" 


n 


n 


s,t^n^ 


( 9 ) 


It is emphasized here that a good point-by-point prediction may 
require better input data than can ever be generated, but the 
point-by-point prediction might be quite informative neverthe- 
less . 


Eq.(8) is the fundamental result for calculating the interior 
sound pressure level. The modal forcing functions are given by 
the term 


Y G (r,H) 



)ip V (x )\p r (x ' )dxdx ' 


( 10 ) 


X X ' 

and must be evaluated for the propeller noise excitation. Once 
a particular harmonic H is chosen, the function TQ(r,H) must be 
computed for all structural modes r. The calculation of ^(^(r/H) 
using the output from the ANOPP program [Ref. 5] is considered in 
Section 3.4. 


High Frequencies 

When the acoustic and structural modal densities are large 
enough, the transmission can be expected to be dominated by 
modal response close to the excitation frequency. Only the 
modes lying in a narrow band Aco containing need be of 
concern. Thus at sufficiently high frequencies the space 
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average mean interior pressure for harmonic H is (assuming 
that only the sidewall contributes) 
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H 

s , t 


i»c 2 
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T ML T t u H 


neAoj n 
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+ (nl + n ' ' ) 
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(ID 


To facilitate the calculation of Eq.(ll), it first is rearranged 
(keeping in mind that 4 , G (r,H) is computed for high frequencies 
for each structural mode as it is for low frequencies, but the 
calculation is limited to those modes whose resonance frequen- 
cies lie close to . Define average loss factors for the 
modes near co^ , i.e., let 


n 


* + n, 


n 


- T1 


n r is the average loss factor for the structural modes. 


Also, let 


s n ( W = [(i-sfj + S n ] 

<MVV ■ + V] 


-1 


Q H " 
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V 2 


(mAV H K 6 
\T " ) T ML T t w H 
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Then Eq.(ll) can be written as 


<P 


i>s,t = ^ E 

neAw n rcAw 


f ' (n s r)g r 'i ; G (r,H) 


2 , , 4 


ITw 
r r 


(12) 


The interior response is sought for the case where w n and w r 
are given equal probabilities of lying anywhere in the band Aw. 


Eq.(12) can be simplified somewhat by noting that the resonance 
frequencies are approximately equal to the harmonic frequency 
w„ of concern. 

n 





W 


H 


£ 

ncAw 


n a n 


E 

reAw 


2 


f’(n>rj gJ (r.H) . 


( 13 ) 


The expected value of the above is now calculated given w n and 
as independent random variables uniformly distributed across 

Aw. 
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In Eq.(l4) } N r is the number of structural modes in Aw and 
is the number of acoustic modes in Aw. 
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Now 


E [g "I = / g ( (!)„/ co ) p(w ) dw 

L n J J n H n ' * n . n 


where p(<o n ) is the probability density function for w n in Aw. 

The probability distribution function is taken to be simply 

0) — w IT 

rw n _ n H 1 

P(w ) = — + o' • 

n A 2 

Aw 


The probability density function is therefore simply 


p<« n > * 


6P( % ) 1 


E [g n ] * ST / %iV“n )<J "r 


that is E[g n ] is just the average 


.Aw 
o> H + 2 


E ^- Sn ^ Aw J Aw 


w H-2 U - V /<) 2 + 


which becomes 


2n n Aw 


H . -1 

tan 


tan' 1 (^±L-) 

VnV VlnV 


Assuming that the product n w„ << Aw (this can be assured by 

n n 

taking the sampling band wide enough), the term in the brackets , 
[ ], above reduces to u , giving 


H - £ 


2n Aw 
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Similarly, it can be shown that 


EtgJ = -Z- 


2r) Aw 
r 


Substituting into Eq.(l4) gives 
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H nsAw 


r ritL K ' )\? n r / J 


But n = N /Aw is the modal density of the acoustic modes, and 
n n 

n = N /Aw is the modal density of the structural modes, 
r r 


Therefore 
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where 
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It can be shown that 
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Substituting this result gives 
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( 16 ) 


In (16), Aco is taken just wide enough to assure smoothness of 
the average calculated in the brackets. 

3-3 Calculation of the Noise Reduction 


In addition to the calculation of transmission of propeller 
noise, there may be occasions where the noise reduction of the 
fuselage is of interest, that is, for a condition where the 
exterior is bathed in a reverberant (diffuse) acoustic field. 

For such a case, the expression for the inflowing power is again 
obtained from Eq.(6). However, in this instance, the joint 
acceptance can be reinserted and brought outside the integral as 
in Eq.(2). Then after the integration is performed, it is found 
that [2] 
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(17) 


where A is the cylinder surface area and A is the coupling area, 
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For n or r = j. 


In. = In 
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(1+c /2 ) 4 to 4 +b . ( 1+c /2 ) 2 to 2 + c . j 
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to J to J 
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where when j = r, r\ . above = p ’ + n ’ ' , and when j = n, n . = n 

J 1 r J ] 

Also 

D nr * (c r- c n) ! + <VV ( Vr ' V„) . 


b n ’ - 2 “n i b r = ‘ 2 “r > 

c n * “n (1 + ; c r ' “r t 1 + <T| r + r >r'> ! ] • 

The above Is Identical to Eq.(5) of Reference [2] except that 
f'(n,r) replaces f'(n,r) and n^+n^!, replaces 

Noise Reduction Calculation 

The interior pressure is obtained by setting W in = W diss (using 
Eqs.(17) and (5)) solving for the mean square modal pressures 
<p 2 > , one-by-one, and adding according to Eq.(3). 

n S y O 

The result is found to be 
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where the exterior space average mean square pressure <p 2 > g t is 
(for the reverberant field) related to the mean square blocked 
pressure by 


<p 2 > , 

F e s , t 


= <P^>/2 


High Frequencies 


When the acoustic modal density is greater than, say, ten or 
more modes in a third octave band, a high frequency model is 
needed, developed along the lines of [1,2, and 3]. In the case 
of an added trim, there are three surface areas of concern: 

1) the transmitting area with trim, which will be called A^ , 

2) the absorbing surface area S, and, 3) the transmitting area 
without trim A'. The total transmitting area is A^+A' and 

the inflowing power takes the form (assuming that only the side- 
wall contributes) 
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where A denotes the surface area of the complete cylinder, 
ff is the sum of the average external radiation loss factor. 


r 

i . e . 


q eX jj and the average structural loss factors q 1 +n " . 
rad r r 


t , is the trim transmission coefficient, and is the resonance 

t n 

transmission coefficient for the diffuse field case, 
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t R = ( 2 V !c o"raaV"d /p o“V) • 1 | P < V A • 

In this expression n^ is the modal density of the fuselage 

structure (modes/rad/sec). The summation r<Am in Eq.(19) 

implies a sum over all structural modes w that are resonant 

r 

below the frequency band of concern. p Q is the exterior air 
density. The radiation loss factors are 
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modes resonant in the band. 


where < j£> r is the joint acceptance averaged over the structural 


Let = field incidence transmission coefficient for mass 
controlled panels as defined in [3]. 
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Then Eq.(19) becomes 
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The power absorbed on the wall is again given by Eq.(5). Since 
at these high frequencies the response is resonant acoustic, co n 
can be considered to lie in Aco and Eq.(5) reduces to 


W„. = jl v 

diss 2 

pc o neAw 


- <P 2 > * • 

s,t 


2k 


U) 



Further, noting that for this case co n - co for any n, and letting 
the group of acoustic modes have an average loss factor h n , this 
reduces to 


W diss ~ P c q <P n > s,t 


w \i "P^s.t 


which is Eq . ( 7 ) , Section III of [2]. Finally, setting 

5c o S 
= ~TuV 


where a is the band average absorption coefficient and S the 
absorbing surface area, the above reduces to the familiar form 
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( 22 ) 


S includes the area A covered with trim and any other absorbing 
surface area. Equating W in to W^j_ ss gives the desired high 
frequency result [^]: 

"Pe's.t “ S + [ T f + T R ( V r, rad ) ] A ' 

2 ) \ J ) 

<p.> , T, T A, + lA T 

*1 s,t t t 

where t = + t r . 


This is the fundamental result with trim present. Note that if 
all transmitting surface is covered with trim. A’ = 0, and 
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On the other hand, if trim does not exist, = 1, A t 

< Pe > s,t = 5S + [vVV' n rad > ] A ' 
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which is the result in [3] without trim. 


Finally, 5 is estimated with the relation 
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W 

c S 
o 



(2H) 


When trim is present n n is calculated with results developed in 
Appendix A. 


3 • ^ Calculation of the Generalized Forces for Propeller 
Noise Excitation 


Consideration is now given to the various terms that appear in 
the equations representing the solutions of the sound trans- 
mission problem. In this particular section, the calculation of 
the generalized (or modal) forces for propeller noise excitation 
is discussed. Specifically, consideration here is given to the 
term ¥ G (r,H) as defined in Eq.(10) and that appears in Eqs.(8) 
and (16). 

ANOPP Computer Program Output 

In the present case, the fluctuating pressure field acting on 
the fuselage due to the rotation of a propeller is determined 
with the NASA Langley computer program ANOPP [5]. This particu- 
lar program is one of a number of such programs to exist in the 
aircraft industry. It was especially developed to allow pre- 
dictions of noise from propellers that have advanced blade 
geometries. In the present case, the concern is more with the 
standard type of general aviation propeller. The ANOPP program 
easily handles the typical general aviation aircraft propeller 
b lade. 

For purposes of the present model, the primary concern is with 
the form of the output data from ANOPP. Referring to Fig. 1, 


26 



on the fuselage structure, there will be a finite number of grid 
points at which the pressure time histories will be available. 

At each point, the modal deflection will also be available, 
say \p (x^). Associated with each grid point is the area AA where 
(Fig.l) AA = A 2 with A being the grid spacing. Now in Eq.(10) 
dx is an element of area, so in discrete form (10) becomes 

V r > H) = i'E E, • 

m m 


where 


* r (x m ) = * r (z k ,8 t ) 




and m identifies each grid point (k,0,k being the axial index 
and H the circumferential one. In the present model, the grid 
is confined to the side of the fuselage where the propeller is 
located. 

To compute 'F p (r,H), the .form of C D , . (x lx ,) that will result 
from the use of the data output by ANOPP is first needed. At 
each point m, there will be a pressure time history that will 
repeat itself over a period T Q = BPF as shown in Figure 1. 
Naturally from point-to-point on the grid, this pressure signa- 
ture will be different. Note that T 0 is 1/B of the period T 
of rotation: T = 60/N, where B is the number of propeller blades. 


Now C_, _(x lx ,) is defined in terms of the pressure time 
Pbl m 1 m 

histories for the point pair (m,m f ). To calculate Cp fel (x [ x , ) 
the equations appearing in Section 3.2 are considered in 
sequence. First 
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R Pbl ( ^mi X m- T) = T- T J /2 Re[p b 1 ( x m ’ } ]Re[p bl ( X m ’ ’ t+T } ]dt (25) 
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I 


Since and p bl ( X m' are P eriodic with period T 0 , 

the time average can be performed over one period, yielding 

T 

r° 

R Pbl (x rJ x m' ;t) = I Re t p bl (x m ,t)]Re[p bl (x rn , ,t+ T )]dt . 

o 

In the present case (x ,t) ls real so 

T 

_ if 0 

R Pbl (x rJ x m’ ;T) * T - p bl (x m’ t>p bl (x m" t+T>dt • (26) 

° J o 

Also since p bl (x m ,t) is composed of discrete frequencies, a 
Fourier series representation can be used 
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Then 


p bl^ x m ,t ' ) 2 + 53 A H cos 


(33) 


p bl ( 
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Substituting these into Eq.(26) (using T q = T^) gives 
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There are four terms appearing in Eq.(35). It is easy to show 
that the second and third terms are identically zero. The 
remaining two terms give upon setting = ^ ~ * 
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It can be seen that the mean. square pressure at x m is 


R Pbl <x mlV 0) ' ‘Pbl'V^’t 
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H=I 

(*S) ! , v + K)’ 
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H=1 


Continuing now, with Rp bl ( x^l x , ; x ) defined, the next step is to 
compute the two-sided cross spectrum using 
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Prom Eqs.(36) and (37), it is found that 
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The one-sided cross spectrum is therefore 
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Finally using c Pbi(*ml x m ,; = C Pbl^ x m! x m ’ ^ ’ lt follows 

that 
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Note that each harmonic corresponds to the tone frequency 

n 

in Eq.(8). To compute Eq.(8), a harmonic is first selected. 
For that particular harmonic, from Eqs.(10) and (40) 
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One of the things to note about this result is that all of the 
information that can be garnered from the propeller noise signa- 
tures is used. However, a disturbing feature of Eq.(4l) is that 
for a big grid, the total number of numerical summations re- 
quired can be extremely large. Since Y G (r,H) must be computed 
for each structural mode, the computations could be very time 
consuming. Fortunately Eq.(4l) can be reduced to a more effi- 
cient result using the trigonometric identity 
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(42) 


In Eq.(42), (2m + 1) summations are required as opposed to the 
m2 summations needed using Eq.(4l). 

Geometrical Considerations and Reflecting Surface Effects 

The ANOPP program computes the acoustic pressure for the case 

of a propeller in a free-field. Here the presence of the fuse- 
lage structure needs to be taken into account since acoustic 
waves will be reflected by the surface. In this model, the am- 
plitudes of the various harmonics are increased in proportion to 
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the incidence angle. The phases computed with the ANOPP 
program are not modified. Until now the amplitudes A™ in Eq.(42) 
have been considered to be the ANOPP results. However, hence- 
forth, they shall represent the blocked pressure amplitudes. To 
distinguish between the blocked amplitudes and the free field 
amplitudes, the free field values from ANOPP will be redefined by 
placing a bar over A™ . The blocked amplitudes are calculated 
from the free field values using 




0.00022^e 0 - 08y ] 


A m 


(43) 


This empirically fitted curve corresponds to that found by 
Magliozzi [6]. Y is the "incidence angle", that is, the angle 
between a line running from the propeller hub to the point x m 
on the fuselage and the normal to the surface at that point. 
The angle is calculated using (refer to Figure 2): 


Y = cos 


-1 


(r cos<f>-acos0 ) ( cos0 ) + (r sin<}>- asin0)(sin0) 

— 2 — — B 

14 rpCOs<}>-acos0 ) 2 + (r^sinf}) -asin0) 2 + (Zp-z) 2 


(44) 


3.5 Interior Coupling Factor f*(n,r) 

This term determines the spatial coupling occurring between a 
fuselage structural mode and a cabin acoustic mode and is 
defined in Section 3.1 as 

f'(n,r) = - N> (x)^ r (x)dx, (45) 

aJ a 

where x represents a point on the transmitting structure. 4> n (x) 
is the eigenvector for the nth cabin mode which is obtained 
using a finite-difference technique as detailed in Appendix C. 
^ r (x) is the structural mode shape and its calculation for the 
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present fuselage model, a cylinder with a structurally integral 
floor, is considered in Appendix D. 


Calculation of f'Cn.r) for Cylinder with Floor 


The modes of the fuselage are given by Eqs. (D.22)-(D.26) in 
Appendix D. For the symmetric modes of the shell, these are 

n* 

^(z,0) = +sin^p £ C Mn (-1)n C0Sn6 06) 

n=0 


and for the antisymmetric modes 
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s r 

C Mn are the generalized coordinates defined by Peterson [13], as 

tabulated for an example case in Table D.l of Appendix D. Also 

as shown in Figure 1, 0 = 0 at the bottom centerline, 0 = n at 

the top centerline. is positive along the outward normal. 

s 


For the floor (plate), the symmetric modes are 
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and the antisymmetric modes are 
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where (Figure 1), x is measured from the center of the floor, 
^ is positive in the upward (inward) vertical direction, and 

c J 

Lp is the width of the floor plate, 

Lp = 2a sin 0 Q 


The acoustic mode eigenvector calculated with the finite differ- 
ence technique takes on discrete values on the bounding surfaces 
of the cabin. Consider the boundary point j = (m,n) of Fig.C-2. 
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From Appendix C, Eq.(C.9), the nth acoustic mode evaluated at any 
point X in the cabin is 


* n ( C) - * ql (S> - oosap. * A ( j ) 


(50) 


On the cabin bounding surfaces, i.e., £ = x the eigenvector takes 
the values 


V x) * V (x) = cosa F~ ♦i ( 9 j ) 


(51) 


where <j>.(0.) is the same as 4>-(j) when j is a boundary point, 
a J a 

To include trim on the sidewall (with a bare cabin floor) the 
transmission coefficient is introduced with the result 


f » (n,r) = f ' (qi,r) = f; I 

Jc\ 

'"i _ 


L p +2a(7r - 0 o ) L 


1 
L 

Vi 
2-rr-e 


qiTZ . Mttz, 
Kin- - C 

L 


cos-V" Sln-^dz 


ar(0)(j>, (6)d6 _ 

o 1 


/ 


o 


L /2 

P ^ p (x)cJ) i (x)dx 

-V 2 

(52) 


A bar is placed over f’(n,r) to indicate inclusion of sidewall 
trim. Eq. (52) is an approximation and not an exact result. It 
should be sufficiently accurate for the present needs. It is 
noted that the 4^(8) are the boundary values of the acoustic 
eigenvector, i.e., the same as the 4> . ( 6 . ) in Eq.(51) and \p (0) 
and ^p(x) are given by the summations in Eqs . (46 )- (49 ) . Since 
4^(0) is known only at discrete points 0 = 0. , j =1, ..., n p , 
the last two integrals in Eq.(52) are approximated numerically. 


Let 0-^ . and 0£j define the points on the circumference of the 

cabin cavity half-way between boundary point (m,n) and the two 

adjacent boundary points. The eigenvector 4>.(0.) is then assumed 

a J 

to apply over the range 8 . to 0„.. For the shell, the second 

1 J 

integral in Eq.(52) becomes 
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that is, if both the acoustic and structural modes are symmetric. 
If both modes are antisymmetric, it is found that the integral 
evaluates to 

n* nfe 
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The integral is zero if one mode is symmetric and the other 
antisymmetric. 
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For two symmetric modes, this becomes 
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and for two antisymmetric modes 
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The integral is zero if one mode is symmetric and the other 
antisymmetric. In this case n^ is the number of boundary points 
at which the acoustic mode has been computed on the floor. 
Finally, let 




For the case considered (shell length L = cavity length L c ), 
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If L / L c , f jyj is given by the relations in Reference [3]. 


3*6 Joint Acceptances for Cylinder with Structurally Integral 
Floor 


The joint acceptance function appearing in Eqs. (16), (18), (19), 
etc., describes the coupling between the excitation field and the 
structure, and is defined by the relation 

j ! («) = 1 (fs , (58) 

r 

where A is the excited structural area, S pt)1 (co) is the blocked 
pressure power spectral density, and S pbl ( x j x’ ; to ) is the blocked 
pressure cross power spectral density. 


It has been general practice, when representing random pressure 
fields with spatially decaying correlation, to use a correlation 
function of the form 

Re(S (x | x * ;oj ) } 

C (x | x 1 ;<o) = 2 . ( 59 ) 

p S (to) 
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It is assumed in the present analysis that this function is 
separable in the longitudinal and transverse directions, i.e.. 


(x | x':o) ) = C x (£, &>)C (£,&») , 

where 

£ = x-x* ; ? = y-y ’ 

The representation for a diffuse (reverberant) excitation field 
is well known and is given by 


c x ( 5>m ) = SLinXMI 

C y (c,<o) = 

^where k is the acoustic wavenumber. 


. ( 60 ) 


Consider the structural mode shapes for a freely supported 
(ideal) cylinder: 

MN 


, r,-> •« Mttz fcos N 6 

* < x) ■ * (2 - e) * sln — (sin N e 


( 61 ) 


The joint acceptances for this simple modal system have long been 
available [17]. Now for the case of the reverberant field, the 
joint acceptance functions for the cylinder must consider both 
sines and cosines of the circumferential wavenumber. One can 

= Ji 


write j£(co) = j*(a)) rev ' 


,* N U) rev ' in the form 


+ isTK-.*) - + (62 > 

where represents the cosine and <j> the sine of the wavenumber. 


The joint acceptance for the axial component of the cylinder 
2 rev 

modes , j M (<o), is given by [17]: 


1 2rev 


(w) = I 1 (M) + I 2 (M) + I 3 (M) 


(63) 
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where 


I 1 ^ M) = 2irMkL 


■{Cin(kL +MTt)-Cin Mit-kL 1} 


i 2 (M) = (Si(kL + Mtt )-Si (Mn-kL )} 


1- ( -1 )^coskL 

I _ (M ) = 

3 (MTT) 2 -(kL ) 2 

Si and Cin are the sine and cosine integrals [18]. 


The joint acceptances for the circumferential modes are [17], 
-lrl. 


Jn ( “ 


fl ■ L sln2^tka(y, -y 2 ) 

^ ^ 2i F a(y ~ -y 2 ) cos 2NT,y l COE SNiy^y^y., 


Jn ( ‘ 


f 1 r 1 sin2Ttka(y,-y 2 ) 

"J J ^ Wy TT y -j sin 2N,y 1 sin 2N„y 2 a yi dy 2 , 

o o \ l 7 

where y ± = y/2ira, y 2 = yV2rra, y and y’ being circumferential 
coordinates of the two correlated points, and k the acoustic 
wavenumber. It follows that 


,rev .rev ,rev 

~ J M ( u> , ip ) + j 2 (u. 


N 


N 


$ ) 


.1 .1 


o o 
This reduces to 


n sin2TTka(y -y ? \ 

2 T ka( yi -y 2 ) cos 2S ’( y l- y 2) dy l d J'2 ' 


rev 

Jj (u) = 


/ l _ 

sin27Tkan /, 

~ 2fkan ( X - n ) co 


S 2irNndr| , 


and 


_ r 1 sln2„kan 
° Jo 2*kan v > 
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where n 


- yi-y2* 


Performing the integration results in 


and 


J j Si[2TT(N+ka)] - Si[2ir(N-ka)]J 


1 

1-cos ( 2 TT ka) 


2 (Ntt) 2 

1- ( — ) 2 
x ' N ' 

) 


.rev 



1-cos ( 2fTka ) 

2 7T 2 (ka ) 2 . 


( 64 ) 
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The cylinder structural model in the present case is much more 

complicated and the cylinder modes are given by Eqs. (46) and 

(47) of Section 3-5* However, the joint acceptances for the 

modes of the cylinder with the structurally integral floor can 

almost be expressed in terms of the joint acceptances for the 

simple cylinder because the mode shapes given in Eqs. (46) and 

(47) are in a form similar to that of Eq.(6l), that is, they are 

given in terms of sine and cosine functions. The joint accept- 

ance for the axial component of the cylinder modes, j^ (w) is 

thus the same as for the simple cylinder and is given by Eq.(63)j 

? rev , 

and the desired joint acceptance j^ (co) is given by 


„rev 

(«) - 


rev 


rev 


M 


(<o) 


N 


(«) 


where for symmetric modes, 

oo 

J^T.) - E(c^) 2 


n=0 


( 66 ) 


(67) 


and for anti-symmetric modes, 


j 


i^> 



( 68 ) 


Note here that N is a counter and is defined only in the sense of 
the duo (r,M), that is, for each r, there is an M and a sequence 

S Y* 

of Cjvjh that are identified with index N. 
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It is fairly straightforward to show that 


j * (“,<l») 


Si(2nka) l-cos( 2u ka) . 

^k^“ ‘ 2Tr 2 TkiT 5 5 

< 

I 0 (n) + I-(n) - I,(n) ; n>0 

2 3 1 


and 


j ^ ( oj , ) = I 1 (n) + I 2 (n) + I^(n) 


( 69 ) 


( 70 ) 


where 


I-, (n) 

I 2 (n) 

I 3 (n) 


1 


4-rrnk ( 2ira ) 


Cin[2ir (n+ka) ]-Cin[ | 2tt (n-ka ) i J 


1 

4nka 


5 (nrr ) 


Si [2 it (n+ka) ] -Si [2 it (n-ka) ] 

l-cos2Ttka 
1- (ka/n) 2 I 


(7D 

(72) 

(73) 


3.7 Resonance Frequencies 


The frequencies « n and &J r of concern pertain to cabin and struc- 
tural modes respectively. The cabin resonance frequencies, co r 
are determined with Eq.(C.7) or (C.ll) of Appendix C using the 
eigenvalues from a two-dimensional finite difference calculation 
as detailed in that appendix. The structural resonance frequen- 
cies, cor> are obtained from the analysis considered in Appendix D 
leading to the associated program output as shown (for an example 
case) in Table D-l . 


3 . 8 Loss Factors 


There are a number of different loss factors that must be 

considered. Referring to Eqs.(8) and (18) of Sections 3*2 and 

3*3, there are the loss factors n n , n r , , and n^' . In 

Eq.(l6), the loss factors n and n appear. Finally, in Eqs.(19) 

i nt 0 x t 

and (23), there are the loss factors n r , ad and h rad . 
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3*8.1 Acoustic Loss Factors n n and h n 
Bare Fuselage (Cabin) 

The acoustic loss factor h n of a bare cabin is defined herein to 
be a measurement for a particular acoustic mode and n n the 
average acoustic loss factor measured for a group of modes 
whose resonance frequencies lie in a narrow band Acd. 


Calculation of n n when sidewall trim is present 


The calculation of the acoustic loss factor at any frequency co 
has been discussed in Appendix A (resulting in Eq.(A.12)). That 
result has to be formatted properly for the case where the floor 
is present. Also, the absorption of the ends of the cabin has to 
be considered. In cases where the sidewall and the surfaces at 
each end of the cabin are covered with a trim, the loss factor 
for the interior is calculated with 


wc e 
o n 


■n 


co ' 


n 


V 


“■ 

“ 

U 4>£dx + 2C e 
J curved j 

£ 4>n dx 

. surface 

ends 


(7*0 


Here E, is the sidewall conductance and F is the conductance of 

e 

the end surfaces. The first integral in Eq.(7*0 is 


/. 


^dx = | 


curved 

surface 


Jib 

e a / \ <j> ? ( 0 . ) 

q fa V j' 


(Vhj] 


(75) 


where n b is the number of boundary points. The second integral 


is 


n. 


f <p*dx = <t>?(j )h 2 C(j ) 

“ends n j = 1 


(76) 


where n^ is the number of interior and boundary points, 
fore at low frequencies, the loss factor is found to be 


There- 



( 77 ) 


'n 


cue 

o 



'n 


V 



z a 

q 


b 

E 

j = l 


»l ( 3 ) [ e 2 J - 9 ij] 


+ 2C e a ! t G d)j 

where $ G (i) is the generalized mass for the two-dimensional mode 
indexed i (i = 0, 1, . .., i’) of Appendix C and z n is the 
normalization constant given by Eq.(C.lO) of that appendix. 

Also e = 2 for q = 0; 1 for q > 0. 

q 


At high frequencies, the average loss factor is used and 

% = TO [ iS + 2 ^e S e] 

S and S e being the sidewall and ends' surface areas. 


is 


(78) 


In cases where the average absorption coefficient a is available 
for the various frequency bands, h R can be calculated with 
Eq . ( 2 4 ) . However, in general it is a that is to be calculated 
from n and ^ is to be analytically derived using Eq.(74) where 
the conductances £ and £ e are based on the trim admittance calcu- 
lations made with Eq.(A.16) of Appendix A. 


3.8.2 Structural Loss Factors 

The structural loss factor n r is defined as the sum of two com- 
ponents, one due to dissipation in the structure (in vacuo) and 
the other due to radiation losses, i.e., 

struc , rad / 70 \ 

n r n r + n r * (79) 

struc 

n r is either a measured value or an estimate based upon 

measurements for, say, similar structures. In the absence of any 
actual measured data, it is recommended that one use [19] 
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struc 


( 80 ) 


where f r is the resonance frequency and the constant C typically 
ranges between 2 and 4. 


2? 

The acoustic loss factor n r is the external radiation loss 
factor for mode r. It is calculated using 

rad 


2p wA rev 

-i£r 


,rev 


o 


(81) 


where j r (w) is given by Eq.(66) of Section 3*6. 


Influence of trim on structural damping 


If trim is present on the sidewall, the loss factor n' r is given 
by Eq . ( A . 7 ) of Appendix A, i.e., 

\h 


|C | 2 
1 w 1 

m 2 oj 4 
r 



(82) 


C w is a (complex) parameter dependent on the values of the 


coefficients of the trim transfer matrix and it is given by 
Eq.(A.8) of Appendix A. 


Influence of Internal radiation: closely coupled structural and 
acoustic modes 


The loss factor modifies the damping level of a structural 
mode when its resonance frequency lies very close to a resonance 
frequency of an acoustic mode. The increased damping is due to 
power flow from the structural mode to that particular highly 
receptive acoustic mode. It is given by [7] 


pc 2 A 2 e f 
K o n 

oj 2 M V 
n r 


* 2 


(n,r) 


[a- 


n 


> 2 , 
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b7(ii 2 ) 
n 


+n 


n 
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(83) 
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The mode r is typically a structural mode lying within a 
bandwidth Aa> n = h n ^ n centered on w . 


Average radiation loss factors 


g x t int 

The radiation loss factors n rad and n rad are given by the 
relations in Section 3-3, namely 


ext 

n rad 


2p <oA 
o 


Time 

o 




rev 

r 


(84) 


int 2pmA ,. 2 / s^rev 

rad Trmc J r r 

o 


(85) 


Again <j-i(^)> T , is the joint acceptance averaged over the 
i i rev 

structural modes resonant in the band. Each j (co) is calculated 
with Eq.(66) of Section 3.6. 


3.9 Validation Studies 


Comparisons of analytical model predictions with measurements 
taken in a model validation experiment are considered in 
Appendix E. 
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APPENDIX A 


SIDEWALL TRIM: TRANSMISSION AND ABSORPTION MODELS 



APPENDIX A - SIDEWALL TRIM 


Transmission 

In Section 3, the basic power flow and dissipation expressions 
used to solve the sound transmission problems are presented. 

The relations include a transmission coefficient t and a 
structural loss factor that appear once fuselage sidewall 
trim is assumed to be present. In this appendix, the analyses 
that lead to these terms are briefly reviewed. Most of this 
material can be found in greater detail in the Phase II report. 
Ref. [7]. As stated in [7], to incorporate the effects of side- 
wall trim in the transmission analysis, the interactions of the 
exterior pressure field, fuselage structure, insulation and 
lining, and interior field must be considered. To integrate the 
trim dynamics into the analysis, the basic expressions that 
describe the response, transmission, and absorption characteris- 
tics of the various components of the sidewall system must be 
considered. To this end, the trim is assumed to be represent- 
able by a transfer matrix: 

a il a i2 c P 1 > (A.l) 

- a 21 a 22-U w l 

where referring to Figure A-l: w-]_ is the displacement of the 

fuselage skin at x, p^ is the pressure on the inner surface of 

the skin at x, w~ is the displacement of the trim panel at x, 
i d - 

and p 2 is the pressure on the trim inside surface at x (cavity 

side). Reference [7] should be consulted for more discussion of 

this trim model. 

In addition to the trim equations there is the basic expression 
for the response of the sidewall. For harmonic excitation it 
is 
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Exterior 


FIGURE A 



-1. A BASIC SIDEWALL TRIM: 
INSULATION AND LINING 
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w 1 ( x ) = /g(x|x';o)[ p 0 (x’) - pj(x')] dx' , (A. 2) 

where G(x|x';co) is the structure's Green's function, and p°(3c') 
is the exciting exterior pressure field. The integral is per- 
formed over the excited structural area. 

Also, the interior acoustic field acting on the trim lining is 
p 2 (x') = - pco 2 jG p (x 1 x' w 2 (x) dx , (A. 3) 

where Gp(x|x';w) is the Green's function for the cavity and p is 
the density in the cabin. 

Finally, the exterior field is 

p°(x') = P bl (x') + P 0 OJ 2 J G°(x| x' ;w)w 1 (x)dx 

= P bl (x') + p r (x' ) , (A. 4) 

P Q is the exterior air density, 

G°(x|x';w) is the exterior space Green's function, 

p r (x') is the radiated pressure field on the exterior surface, 

Pp-r(x') is the blocked pressure field on the exterior surface. 

Equations (A. 2), (A. 3) and (A. 4) and the trim transfer matrix 
(A.l) form a system of five equations in five unknowns: p^, p 2 , 
W]_, W 2 , and p r (or p°) ; the blocked pressure field is assumed to 
be determinable once the geometry of the structure is fixed and 
the excitation defined. The solution of this system of equa- 
tions determines the effects introduced by the presence of the 
trim. Note that if the trim is removed pij- = p^ and w_^ = w^ 
reducing Eqs.(A.2), (A. 3), and (A. 4) to the set solved in 
References [1] and [2]. 
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In Reference [7] it is shown that the trim panel displacement is 
very nearly given by the result 


w 2 (x) = c'^Cx) 


(A. 5) 


where C w = a ^ and a^ comes from- the trim transfer matrix. 


The 


acoustic power radiated into the cabin is computed with the 
theory presented in Ref.[l] and [2], and it is shown in Ref. [7] 
that the inflowing power is given by Eq.(6) of Section 3*1 of 
this report, where C w appears in the form of its magnitude 
squared,- which is defined to be 


= C 


w 


2 = |a 


, - 2 


11 


( A . 6) 


Since this term multiplies a power flow expression that basical- 
ly remains the same whether trim is present or not, I C w I is 


easily interpreted as a transmission coefficient for the trim. 


It is also found that when trim is present, the structural loss 
factor n r is augmented. Thus if n r is defined to be that due to 
dissipation in the structure when trim is absent and is the 

total structural loss factor when trim is present, then accord- 


ing to the analysis in 

C 7 ] 

j 
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2 2C 1 

n t 2 _ J 

w 
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2 4 

m oo 
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-E + n 2 

o > 2 
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(A. 7) 


where 


C 

w 


'w 


+ i C 


w 


i 


and is obtained from the remaining three elements of the trim 
transfer matrix from the result 



Trim Transfer Matrix 


The transfer matrix for a layer of insulation is the following 
(Appendix A of Ref. [7]): 


C\J 

a 


"C -ws ■ 


,p l 



s „ 



t V 2 


w c 


V 1 , 


where C = cosh YL ; S = sinh YL ; Y = a-ik = a -i2ir/Am 
Y is the propagation constant of the insulation (complex). W is 
the wave impedance of the insulation (refer to Fig. A. 2). It is 
noted here that Eq.(A.9) is not in the form of Eq.(A.l). The 
matrix of concern, given in Eq.(A.l) relates pressures and 
displacements, rather than pressures and velocities. Since 


v i = -iww i 

Eq.(A.9) is more properly written as 



'C + lows' 


f P l| 


-IS 




. coW . 


l w lj 


(A. 10) 


In the case where the trim is nothing more than a layer of in- 
sulation, the trim transfer matrix would be given by Eq.(A.lO). 


If there is a lining consisting of a limp mass, the transfer 
matrix across it is the following: 



'X ai 2 m t +ia. I n t m t -| fp 


w n 


(A. 11) 


where the index 2 is the output terminal of the four-pole and 1 
is the input. m t is the mass per unit of area of the lining and 
n t is some assumed loss factor that arises because of flexure of 
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FIGURE A-2. ACOUSTICAL PROPERTIES OF OWENS-CORN INC 
PF-105 FIBERCLAS: DENSITY 9.6 kg/m 3 
(0.6 lb/ft 3 ) , [21]. 
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the lining (the stiffness being ignored). For a multielement 
trim the transfer matrix across all elements of the trim is 


fH. 

a ll a l2 


' a ll a 12 _ 


a ll a l2 


a ll a i2 " 

f 1 

l w 2 i 

_ a 2l a 22. 

n 

- a 2l a 22_ 

• • • 

n-1 

- a 21 a 22. 

2 

- a 2l a 22_ 

l w l 

1 


a ll a l2 


a 2l a 22 




i w i 


(A. 12) 


where l is trim element in contact with the skin and n is the 
finishing element in the cabin. The transfer matrix for the 
trim of Figure C-l can be determined using Eqs. (A. 10), (A. 11), 
and (A. 12) 

(A. 13) 


’P' 


iW, 


( C-icom^^+wri^m^y-) (iwWS+co 2 m^C + ia)' i n^m + C ) 


' t t 


IS 

"cjW 


According to the definition of the trim transmission coeffi- 
cient, see Eq.(A.6), for this case 

, -2 


T t = 


111 


where 


a ll ' (C + cor, t m tW ) " 1Wm tw 


(A. 14) 


Note that C, S, and W are complex. 

Eq. (A. 12) can be used to as high a frequency as is found 
(through experimental comparison) to be valid. In practi- 
cal trims, the mechanical vibration transmission from fuselage 
skin to trim panel will become significant at high frequencies 
and eventually the errors involved with (A.l) may require that 
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some alternate approach be considered to obtain the correct 
input data for the trim. Because of the result in Eq.(A.5) one 
can always take the ratio of say the measured mean square dis- 
placements (or velocities or accelerations) of the fuselage skin 
and the trim panel, average the results over the transmitting 
surface and estimate ]C W | 2 (and thus t^). Also by measuring the 
actual damping of the fuselage wall with trim installed, can 
be determined and basically all of the pertinent information 
required from the terms of the trim transfer matrix is obtained. 
This.allows a fairly elementary approach to the use of any 
available experimental data or of any analyses that have focused 
on the vibration transmission (or shorting) problem. 

Sound Absorption 

In addition to the transmission and damping effects, the trim 
installation will be a sound absorbing system, that is, it will 
take energy from the cabin space. It is shown in [7] that this 
absorption capability is also describable from the terms of the 
trim transfer matrix. In fact, the loss factors for the cabin 
acoustic modes can be determined with the relation 

ojc e f 

h n = IT 2 / 5(x)<f>£(x)dx , (A. 15). 

x 

where ?(x) is the conductance looking from the cabin into the 
sidewall treatment. For a case where K is independent of x, we 
can use 


K = Re[8] 


where 8 is the admittance given by [7] 


8 = pc 


-0j2a 21^1 + ■*- a>a 22 
lcoa ll Z l + a 12 
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(A. 16) 



It is seen that the transfer matrix (A.l) describes the trim 
absorption, except for the presence of the term which is the 
structural impedance. In [7], a simple model 


Z 1 = p o c o“ iam 


(A. 17) 


is recommended, and this is used in the present program. 


Vibration Transmission to Trim 


As stated previously, the mechanical vibration from fuselage 
skin to trim panel can become significant at high frequencies. 
According to the results found in Ref. [7] and summarized by Eqs. 
(A. 5) and (A. 6), the trim transmission coefficient is simply 

' 2 

<wi> 

t = [C w | 2 = — — > (A. 18) 

<w 2 > 

that is, it is the ratio of the mean square trim panel displace- 
ment upon the mean square skin displacement. As long as trans- 
mission through the Insulation is dominant, = la 11 l ; but 
once the mechanical vibration transmission is dominant, <w 2 > 
will be larger than that estimated with Eq.(A.6) and instead, 
<W 2 >/<w^> must be a measurement or an analytical prediction for 
the mechanical vibration transmission. At high frequencies, a 
simple expression for this ratio can be derived using statisti- 
cal energy analysis procedures. The transmission coefficient 
for the mechanical path is [20]: 


m 

T t 


m t A t n t 

mAn 


n 21 

11 21 + ’H 


(A. 19) 


where m^_ is the trim lining mass per unit area, A t is the trim 
area and n, the trim modal density, m. A, and n are the same 

T/ L 

quantities for the fuselage sidewall. n t is the trim loss factor 
(in flexure) and is a coupling loss factor for transmission 

from trim to skin. 
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For N point connections between skin (frames) and trim [20]: 


' h t c t\ (mhc i )(m t h t c t ) 


F) 


(A. 20) 


co > (mhc^ + m^h^c ^.) 2 


where h and c are the trim lining's mass per unit area and 
t t 

dilatational (longitudinal) wavespeed; *h and c^ are the 
corresponding skin properties. For a line support, the coupling 
loss factor per unit length is [20] 


• ( 4 )‘ h ( 4 *-) 




(mh^c|* + m^h^c ^) 2 


(A. 21) 


It is emphasized here that the finding in Ref. [7], namely 
Eq.(A.5), forms the theoretical basis and justification for this 
approach to trim transmission. Using these results, it can be 
stated that for any trim, the i ^ to be used in Eqs. (8), (16), ( 1 8 ) 
and (23) is given by 


*n r 2 + 


(A. 22) 
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APPENDIX B 


FUNDAMENTALS OF PROPELLER NOISE THEORY 



APPENDIX B - PROPELLER NOISE 


The goal here is to develop the basic equation used for 
propeller noise prediction in the ANOPP program. There is no 
attempt to discuss its mechanization, although appropriate 
references for further study are included. Consideration is 
limited to propellers with subsonic tip speeds and to the thick- 
ness and lifting noise components of the sound generation. 

As stated earlier, there is always a problem determining just 
what level to begin a technical presentation when there are a 
number of different disciplines involved. In this particular 
case, since propeller noise is not an easy subject to under- 
stand, the authors have taken what might be considered the path 
of least resistance by utilizing a result by Goldstein that has 
been described as the fundamental equation of sound generation 
in the presence of solid boundaries [Eq.(3.6) of Ref. 8], namely 


(? 


• iff 


9 2 G 


ij 
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In Eq.(B.l), p'(r,t) is the fluctuating fluid density at the 

observer location r at time t. The equation applies to any 

region V (t Q ) which is bounded by surface S 0 (t Q ) in arbitrary 

motion. S (t ) can be an inner bound as in the case of a 
o o 

propeller blade. There are three terms in (B.l). The first 
involves Lighthill's stress tensor T. . . This term ultimately 
leads to sound of the quadrupole type being radiation emitted 
from the volume of fluid around the blade and is important only 
when the tip speeds are in the transonic region. The remaining 
two terms are the lifting and displacement noise components 
respectively. P-Po is the difference between the pressure on 
the blade and the ambient pressure, n^ is the ith component of 

g 

the unit normal to the blade surface, is the ith component of 
the blade surface velocity, and G is the free space Green's 
function. 


G(? o 


t Q ) |r,t) 


1 

Wr 


S(t -t + R/c ) 


(B . 2 ) 


Here r is the source location and R is the distance between 
source and receiver locations 


R = i r ~ r 0 i = |x-y| 

Note that the normal gradient is 

3 • 9 

3n = i 3y 1 

o i 

where y i is the ith source coordinate and a sum over the index i 
is implied. The normal surface velocity is 


V 

n 


Tr s 
n. V. 
i i 


~ ?ts 
n*V 
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where n is the unit normal vector. 
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One reason for beginning with this result, is that it is de- 
veloped carefully, unambiguously, and ingeniously by Goldstein 
[8]. A reader can easily wade through the algebra that leads to 
Eq.(B.l). The remainder of this appendix is concerned with the 
manipulation of this result to obtain the ANOPP equation, what 
shall be called here as the propeller noise formula. 


To begin, it is noted that the relation between the fluid den- 
sity fluctuation and the fluid pressure fluctuation, (neglecting 
heat conduction in the fluid) is simply 

P* = c*p ' 


For convenience, the prime on the p is dropped on the left hand 
side, and on the right hand side, we replace the airfoil pres- 
sure (p-p Q ) by l. Now p is understood to be the fluctuating 
pressure and it follows that Eq.(B.l) can be rewritten as 
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(B. 3) 


In the above R is independent of time. 
Consider first the evaluation of 
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Now 


where 


Also 


where 


and 
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u r is the so-called unit radiation vector 
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latter result follows because R is independent of t 0 . 
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Following the general approach used in propeller noise predic- 
tion, it is now convenient to introduce a Lagrangian coordinate 
system (so-called n system [9,10]) that rides in the body and in 
which S G (t 0 ) remains fixed (see Figure B-l). 

The velocity V and acceleration a of any point n in this coordi- 
nate system is 

^ ~ Ft ^ #t o^fi = fixed 
o 1 

a = ( n , t 

9t * o n = fixed 

Since S 0 (t 0 ) = S Q is fixed in this system, i.e. , S Q is a rigid 
surface; every point has velocity 

V = V 0 (t 0 ) + S(t o )xn 


where V Q is the velocity of the origin of the n system and co is 
the angular velocity of S Q (or the n system). Since n and y are 
Cartesian, the Jacobian of this transformation, i.e., 

( r Q or y). fi 

is unity. 


In the f) system, the limits of integration of the surface inte- 
grals are independent of t and the order of integration in 
Eq.(B.3) can be reversed to yield 

f f 1 n*u l 
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Receiver 



FIGURE B-1. LAGRANGIAN COORDINATE (n) SYSTEM 



where now R is a function of time, i.e.. 


y=r o =f+n=n+ E(t Q ) = ? 0 (t Q ) , 

arid 

R = |r-r o l = |x - n - C(t Q )| • 


To evaluate the integrals in (B.^) the following identities from 
generalized function theory are used [9] 



f (t )<5[g(x )] dx = 
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where x^; is the ith root of g(x^) = 0 . 
0 0 


Also 
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The first identity (B.5) Is used on the first term in (B.4) to 
obtain 
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Therefore , 
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t =t e 
o o 


where t e = t-R/c 0 . . 

The second identity (B.6) is used on the second term in (B.4) 
and then the first identity (B.5) is applied to- further reduce 
the second term. Using 


f(t c ) * | [K • V + f> 0 c o fi 
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Now using (B.5), setting 
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gives the above as 


1 3 

| i-M | at 
1 r 1 o 

o ' ”o 

Substituting (B.7) and (B.8) into (B.4) gives, after changing 
the sign so that the normal points out of the blade. 
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Eq.(B.9) can be written in an alternative form which is faster 
to execute on a computer and more accurate [12]. First, the 
acoustic pressure is separated into loading, p^(r,t) and thick- 
ness Prp(r,t) contributions, such that 


p(r,t) = p L ( r , t ) + p T (r,t) 


Then the two components can be written in the forms 
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and 
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Mach number terms and M 2 are given by 


M. = V./c and M 2 = V.V,/c 2 
110 1 1 o 


A /S 


where (u *n). Furthermore, the dot on and £ 

o 


denotes the rate of variation of these parameters with respect 
to time t . 


Equations (B.10) and (B.ll) have been coded at NASA Langley for 
use in ANOPP. The reader should consult references such as [11] 
and [12] for details concerning the mechanization of these 
results. Also, an alternative derivation of the present results 
is provided by Farassat in References [10] and [12]. 
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APPENDIX C - CABIN ACOUSTIC MODES 


The acoustic modal characteristics of the cabin are defined by 
the mode shapes <J> (£), Z being an interior point, resonance fre- 
quencies a) , and the mode normalization given by the constant 
e . In the present case, consideration must be given to the 
determination of these quantities for the case of a cylinder 
with a floor partition defined by the angle 6 q as shown in 
Figure C-l. In the ideal complete cylinder the modal properties 
can be determined in closed form by an analytical solution of 
the wave equation, subject to the appropriate boundary condi- 
tions. This is possible because the wave equation is separable 
in cylindrical coordinates and the boundary conditions can be 
expressed in these coordinates. When the floor is present, it 
is no longer possible to derive the mode shapes analytically 
since the boundary conditions are irregular. Thus it is neces- 
sary to resort to numerical methods. There are two possibili- 
ties, either finite differences or finite elements. Since in 
the present case, the modal characteristics in the axial direc- 
tion are known, a two-dimensional problem remains, and the 
finite difference technique, which is the simpler of the two, is 
chosen. First, the two-dimensional problem is solved. Then the 
axial modal information is factored in. Next the normalization 
of the data is defined. 


Finite Difference in Two Dimensions 


In the cavity (cabin), the Helmholtz equation applies. In the 
two-dimensional problem, using central differences, it is found 
that for the grid of Figure C-2, the pressure P obeys the relation 

m,n m+l,n m-1 ,n” r m,n+l~ m,n-l ~ . r m,n > (C.l) 

where h is the grid spacing, and k =w/c Q . 


C-l 



Cabin: 



FIGURE C-1. CABIN ACOUSTIC SPACE 
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(m,n) = (NX, NY) 


FIGURE C-2. FINITE DIFFERENCE GRID NOMENCLATURE 
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Boundary Conditions 


The boundary condition is that the outward normal gradient is 
zero (the wall admittance 6 is assumed to be sufficiently small 
to allow this assumption). Thus, referring to Figure C-2, in 
the region © Q < 9 <2^— © Q : 


( C . 2 ) 


If -e o <e<e o , the boundary condition is obtained from Eq.(C.2) by 
simply setting 6 mn = o whenever 0 mn < 0 O • 

Solution 


Let the column vector Pj = ^Pj^njxl be tbe pressure at interior 

and boundary points. Also let = {P F ) be the pressure at 

^ n e xl 

the. exterior points, adjacent to boundary points only. Using 
the recurrence relationship for all interior and boundary 
points, i.e. , Eq.(C.l), gives (in matrix notation, with 
X = k 2 h 2 ). 


R i p x + 


r e p e ap i 


( C . 3 ) 


where Rj has dimensions nj xnj and Rg has dimensions njxn e , 
where nj is the number of interior and boundary points and n g is 
the number of exterior points. The matrices Rj and R £ will be 
different for symmetric and antisymmetric modes. 

The boundary conditions (Eq.(C.2)) give (n e -2) equations. The 
tangential gradient is assumed zero at two boundary points close 
to 0 = t r/i| and 3^/4, giving n g equations. When combined, these 
take the form 


Vi + Ve ' 


C-H 


0 


5 


(C.4) 



where Bj has dimensions n e xnj and B e has dimensions n e xn e . 

Since the matrix Bg is non-singular 

p e * - b e 1b i p i : (c -5) 

Substituting this gives 



The eigenvalues and eigenvectors are calculated for symmetric 
and antisymmetric modes, separately. The modes are then com- 
bined, and ranked in ascending order of frequency and the first 
20 modes only are used. The vector P^. corresponding to X is 
the mode shape where i is the 2-D mode counter, and j 

defines the position in the fuselage cross section. The nodes, 
n, , represent only half the cylinder (i.e., X > 0). The values 

u 

of the eigenvectors (on the boundary) will differ for symmetric 
and antisymmetric modes for X < 0. For symmetric modes, 

4> ± ( J ) = ^(X.Y), 


d 1 (X,Y) = <*>.(-X,Y) 

For antisymmetric modes, 


<j) 1 (X,Y) = -<f>.(-X,Y) . 

Each mode must therefore be identified as symmetric or anti- 
symmetric when the boundary values of the’ eigenvectors are used. 
Apart from the boundary values of the eigenvectors, there is no 

need to distinguish between symmetric and antisymmetric modes in 
✓ 

the cylinder response program. In addition, it should be noted 
that the first symmetric mode represents uniform translation with 
zero frequency. 


C-5 



are used to calculate the resonance 


The eigenvalues, ^ , 
frequencies for the (three dimensional) cabin with 


f 


n 




9 


(C . 7 ) 


where h is the grid spacing ( =Ga), and L and a are the cylinder 
length and radius respectively. The frequencies output by the 
program are 



i.e., the frequencies for q = 0 modes and. a = 1 meter. 
Normalization 

The acoustic modes <f> n for the three dimensional cabin are 

V^ )= V X ’ Y ’ z) = V (X,Y,z)=C ° S ^V X ’ Y)=cosa ^i (J) ’ (C.9) 

where ^(XjY) is the eigenvector for the two-dimensional mode 
ranked i (i = 0,1,2,...) as calculated with the finite differ- 
ence technique; i.e., <p^(X,Y) is a finite dimension column 
vector, which contains the values of at all coordinate 
positions j within and on the boundary of the cabin space. The 
normalization of the modes is arbitrary. The maximum value 
achieved at any coordinate positon has been chosen to be unity 
and the other values adjusted to retain the computed ratios from 
point-to-point. The normalization is carried into the trans- 
mission prediction with the parameter -a . In the case where the 
floor is present the integral required to determine -e is 

- | e q £>’(J)h ! C(J) 
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where j counts over all interior locations. C(j) = 1 for 
boundary points and centerline locations, C(j) = 2 for interior 
points. Also = 2 for q = 0; 1 for q > 0. The generalized 
mass for mode i as defined in the finite difference program with 
unit radius is 


*G (1) = 



) G 2 C ( j ) 


(C. 10) 


The volume 
Thus 


enclosed is 

V = a 2 L 1 7T-0 + cose sine 


e = e . = 
n qx 


[■- 


o 


] 


2 | it — 0 + cos6 sin6 

_o o o 

V 1)£ q 


] 


Sample Results 


(C.ll). 


Figure C-3 shows examples of the finite difference calculation 
for a case where the floor 6 q is 56.6°. The first twenty modes are 
shown and are ranked according to the A^. The resonance fre- 
quencies are for a = 1 meter. The 1 meter results given in 
Figure C.3 are simply the f^ in Eq.(C.8). The resonance fre- 
quencies for arbitrary radius a (in meters) is obtained from 
Eq.(C.7)» Eq.(C.7) can also be written in the form 



where fj_ is the 1 meter result. Note that for q = 0, 

f. 

f = — 

Oi a 

Table C-l gives the computed generalized masses, $^(i) for the 
1 meter (unit) radius. 
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MODE 0. SYMMETRIC 
FREQ. = 0.00 


MODE 1. ANTISYMMETRIC 
FREQ. = 99.5 


MODE 2. SYMMETRIC 
FREQ. = 120.4 
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FIGURE C-3. ACOUSTIC MODES FOR A CASE WHERE 
= 56. 6°, q = 0, a = 1 rrieter 



MODE 3. ANTISYMMETRIC 
FREQ. = 175.9 


MODE 4. SYMMETRIC 
FREQ. = 187.9 


MODE 5. SYMMETRIC 
FREQ. = 220.4 
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FIGURE C-3 
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MODE 6. ANTISYMMETRIC 
FRE. = 253.2 


MODE 7. SYMMETRIC 
FREQ. = 259.2 


MODE 8. ANTISYMMETRIC 
FREQ.= 289.7 



FIGURE C-3. CONTINUED 



MODE 9. SYMMETRIC 
FREQ. = 300.2 


MODE 10. ANTISYMMETRIC 
FREQ. = 324.6 


MODE 11. SYMMETRIC 
FREQ. = 343.7 



FIGURE C-3. CONTINUED 



MODE 12. ANTISYMMETRIC 
FREQ. = 356.9 



MODE 13. SYMMETRIC MODE 14. SYMMETRIC 

FREQ. = 366.1 FREQ. = 381.0 



FIGURE C-3. CONTINUED 
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MODE 15. ANTISYMMETRIC 
FREQ. = 393.3 


MODE 16. SYMMETRIC 
FREQ. = 418.4 


MODE 17. ANTISYMMETRIC 
FREQ. = 424.0 
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FIGURE C-3. CONTINUED 



til- 


MODE 18. SYMMETRIC 
FREQ. = 434.1 


MODE 19. ANTISYMMETRIC 
FREQ. = 447.9 


o 




FIGURE C-3. CONCLUDED 



TABLE C-l. EXAMPLE PROGRAM OUTPUT 

Phase III Test Article 

Angle subtended by floor edge with 
vertical, 8 q = 56.6 degrees 


MODE FREQUENCY 


GEN. HASS 


0 

0.0 

1 

99.501 

2 

120.379 

3 

175.389 

9 

137.927 

5 

220.381 

6 

253.162 

7 

259.176 

8 

289.679 

9 

300.189 

10 

329.597 

11 

393.723 

12 

35b. 999 

13 

366.051 

19 

381.030 

15 

393.310 

16 

913.379 

17 

929.029 

18 

939.105 

19 

997.637 


SY/Irt 2.59631 
ANTI 1.03531 
SYrt* .72535 
ANTI .68099 
SYM8 .72900 
SY.Irt .99315 
ANTI .59535 
SYrtfl .50751 
ANTI .51521 
SYMrt .36571 
ANTI .92606 
SYrtrt .95990 
ANT I . 9600 5 
SYM* .36319 
SYn /1 .39391 
ANTI .23999 
SYrt* .39305 
ANTI .37055 
SY*,“! .90299 
ANTI .95990 
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APPENDIX D 


FUSELAGE STRUCTURAL MODEL: 
CYLINDER WITH INTEGRAL FLOOR 



APPENDIX D - FUSELAGE STRUCTURAL MODEL 


The fuselage modal characteristics are defined by the eigen- 

V — 

vectors (x), resonance frequencies and the mode normaliza- 

tion given by the constant M r (the modal mass). In the present 
model, these quantities are computed for a ring-stringer 
stiffened cylinder with a floor partition that is structurally 
an integral part of the fuselage. 

The basis of the structural model is an analysis of the free 
vibrations of a circular cylindrical shell with a longitudinal 
interior plate (Fig.D-1) by Peterson and Boyd [13, 1A]. A review 
of their analytical approach is presented in this appendix and 
then the generalization of their results to enable calculations 
for the ring-stringer stiffened shell is considered. 

Displacement Functions 


The displacement functions are assumed to be finite series. For 
the shell: 


M* n* 


u = 


Z Z U® x M (xH (6) 


M=0 n=0 
M* n* 

■ E E 


Mn uM un 


VlX H (x)i (9) 


mTo n=0 Mn V M vn 


(D.l) 


w 


M* n* 

■ E E 

M=0 n=0 


W Mn X wM (x) «Vi (0) • 


The plate functions u p , v p , w p are similarly expanded. 


The longitudinal functions X(x) are the same for the plate and 
the shell. These functions are expressed in terms of a single 
function (^(x), in the following manner: 
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FIGURE D- 



. CIRCULAR CYLINDRICAL SHELL WITH A 
LONGITUDINAL PARTITION 


D-2 


X uM = <W X ) 


X VM ' V x) (D.2) 

X wM = V x) * 

The functions ^(x) are the mode shapes of a uniform beam. 

The boundary conditions used in this model are those of a 
supported beam, so <f>^(x) = sinMux/L. Note here that the shell 
axial coordinate is x. In the cylindrical coordinate system 
used in the body of this report, the axial coordinate is z. No 
confusion should result from this temporary change in nomencla- 
ture. 

The circumferential functions for the symmetric modes of the 
shell are 

iii = cos n0 
un 

^vn = sln n0 
^wn = cos n0 

For the antisymmetric modes, (D.3) 

ib = sin n0 . 
un 

^ = -cos n9 

r vn 

\b = sin n0 
wn 

A similar set is used for the symmetric and antisymmetric modes 
of the plate (floor). 

Constraint Equations 

The floor partition can be taken to be fixed or pinned, (hinged) 
along the line of attachment to the shell. For a rigid 
attachment, the shell and plate displacements obey the following 
relations 
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(D.H) 



w 

sin0, + 

v cos0. 

= v 

s 

1 

s ■ 1 

p 

w 

COS0, - 

v sinO, 

= w . 

s 

1 

s . 1 

p 

v 

s 

, 3 w 
1 s 

3 w 

+ _ _E = 

0 

r 

r 3 6 

3 y 


In the case of a hinged connection the last equation is dropped. 
Equations of Motion 


Peterson [13] used Hamilton's principle to derive the equations 
of motion. He assumed that the cylinder was of uniform thick- 
ness t and allowed for different moduli! E x and Eg, but no 
stiffeners. Following [13], the displacement functions are 
written in the matrix form 


“ s | 

V s| ' ^ {q s } 5 (D. 5 ) - 

w I 

S / 

where [N„] is a matrix of size 3 x 3M*n* and {q} is a vector o; 

o 

the generalized coordinates of the shell given by 




where 


(u } = 
s 


U 

U 


00 

01 


U * 
y on* 

U ln* 


U M*n* 


h { v s } 


= < 


V 

V 


00 

01 


V 

V 


on* 

In* 


v. V M*n*) 


W. 

W' 


00 

01 




w r 

w 


on 1 


In* 


W 

^ M*n k J 


(D.7) 
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The displacements of the plate are expressed in the matrix form 


= [N p ] {q p } , (D.8) 

where the elements of {q } are the generalized coordinates of 
the plate, arranged in the same manner as the elements of (q s ). 
For either the shell or plate the strain energy can be expressed 
as 

u = \ I (a} T (e}dS , (D. 9) 

J s 

where (a) is a vector of stress resultants and (e) is a vector 
of strains and curvatures for the plate or shell. 


P 

v 

P 

w 


For a linearly elastic material, 

(a) = [D] (e) 


(D.10) 


where [D] is a matrix of elastic constants. The strain- 
displacement relations, obtained from Love's shell theory and 
classical plate theory, have the form 


(e) = [G] (u), where (u) = Jv • 

\w , 


( D . 1 1 ) 


Substituting Eqs.(D.lO) and (D.ll) into (D.9), and using 
Eq . (D. 5 ) or (D.8) leads to 

U = \ £ ([G] [NJ (q)) T [D] [G] [N] { q) dS . (D.12) 

o 

For either the shell or the plate the kinetic energy is (in 
terms of the generalized coordinates) 

U = | {q} T [K](q} , (D.13) 

where 

[K] = / [[G][N]] T [D][G][N] dS . (D.Hl) 

s 

D-5 



Writing the strain energies for the shell and plate as 


°s = K> T[ V<Q s > ana u p =■ |{q p ) T [K p ]{q p ) , 


gives the total strain energy: 


U = U + U =4 
s p 2 


(q_ 1 

T 

1 

O 

J 


fq 

CO 

) 


s 


Is 

q 


0 

K 


|q 

Ip 


L PJ 


Ip 1 


= \ {q} T [K*]{q} 


The kinetic energy of the plate or shell is 

T = 5- / dS. 

d S 

In terms of the generalized coordinates 

T = \ J s p{q} T [N] T [N]{q} dS . 
The mass matrix is defined as 

M = / p [N] T [N] dS , 

S 

giving for either shell or plate 


(D. 15 ) 


T = \ {q} T [M]{q} 


The total kinetic energy for the system is the sum of the 
kinetic energies of the shell and plate, T = T s + T p , yielding 


fq i T rM 

0 1 

rq 1 

J.s U s 


1 . s \ 

lq J 1 0 

M 

lq J 

p L 

P J 

4 P 


[} 1 [M* ] ( q } 


where 


(q) 


■ Cl 


is the combined vector of all the generalized coordinates of the 
shell and plate. 


The components of the coordinate vector {q} are not independent 
because the constraint equations must be introduced to insure 
displacement compatibility at the interface between the plate 
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and shell. The equations expressing this compatibility can be 
written in terms of the generalized coordinates in the following 
matrix form: 

[ G ] { q } = {0} . (D.17) 

Applying Hamilton's Principle and adjoining the constraint 
equations by introducing a vector of Lagrange multipliers in the 
standard way leads to the equations of motion and constraint for 
the system. 

If { q) is partitioned into a set of independent coordinates (q q ) 
and dependent coordinates {q 2 ) 3 the constraint equations can be 
manipulated and the dependent coordinates and Lagrange multi- 
pliers can be algebraically eliminated. The equations of motion 
then take the form: 

[M] { q x } + [ K ] { q x } = {0} , (D.18) 

in which 

[M] = [E] T [M*][E], [K] = [E] T [K*][E] , (D.19) 


where 


[E] 


I 

[c 2 ] _ 1 [c 1 i 


(D.20) 


and Cj and C2 are obtained from the partitioned constraint 
equations 

[ c i: c 2 ]{q*} - {0) ■ < D - 21 > 

The eigenvalue problem is obtained by letting {q^} vary harmoni- 
cally with time to yield finally 

[K]{ qi ) = o) 2 [M]{ qi } . ( D . 2 2 ) 
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Now, the primary concern in this appendix is to modify the Peterson 
results to enable calculations for a stiffened shell. From Eqs. 

(D . 10 ) and (D.ll), 

(c } = [D] { e } = [D][G]{uJ 
s s 

Peterson's results for the above (Appendix A of Ref.[13]) can 
be written in the modified form 


x 


N 


y 


N 


xy 


> =• 


M 


M 


M 

v Xy J 


k 3 
xs 

3x 


v 0 k xfx 


1 k 


. 3 
x 030 

- v k 
r x 0 

3 

0R 3 0 

i k 

r K 0R 


in f— r t 
r xeie x0 3 x 

1 3 * 

0 V D q - 2 7~X~ 
x 0r 2 30 


0 D — -2— 

0 r 2 30 


tj_ 1 3_ 

! x 0 12 r 3x 


n il n 1 3 2 

* xs 3x 2 " x 0 r 2 3 0 2 

v D ii_ i n ll_ 
0 x 3 x 2 " r 2 0R 30 2 


-G • — 

^ __ TO 


2 3 


x0 12 r 3x30 


r ^ 


u 


V 


w 

s 

V. J 


r (D.23) 


In Eq.(D.23), Peterson's k , k D , D , and D Q have been replaced at 

x t) x b 

various points in the matrix by the augmented terms 


xs 


k x + i d / stringer 


k 


0R 


k 0 + 


D 


xs 


D 


0R 


D 0 + 


(“) 

(¥) 

(“) 


ring 


D x' + V d / stringer 


ring 


(D . 2^4 ) 


The result may be compared to Eqs. (25)— (27) on Page 10 of Mikulas 
and McElman [15]. If the assumption is made that the stiffeners 
used by Mikulas and McElman are symmetric about the skin (i.e., 
z g and z r are zero), the equations are very similar except that 
the terms with the asterisks do not appear. 
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The strain energy term is obtained from Eq.(D.12). The result 

T 

of concern is [G] [D][G] and is given by 


, 3 2 ■ 1 r t 9 2 

"xs 9x 2 r r xe 30 2 



CD 

X 

X 

? 

( 1 1 1 

3 2 . 1 p t 3 2 

3x36 r x0 3x39 


h v x k e W 


L 

- 


- _ 




1 3 2 

+ - G n t 3 


r 0 x 3x90 r xe 3x30 


k. 


^ + G _ A t 


r x ^0R 30 2 X 0 3x 2 


+ D i_ ii 2 _ + 0 1 J)jL_ 

6 r 4 98 2 x0 12 3x 2 


r 2k 0R 30 “ r 2 V 0 D x3x 2 30 


.In t-G^Ll 3 


r 4 8R 30 3 ” x0 12 r 2 3x 2 30 


r V 0 k x 3x 


— k 3 v o 


r 0R 30 x 0 r 2 303 x 2 


1_ n 3 3 2_ r tii 2 

r 4 u 0 SfF r 2 x9 12 Sx^M 


k„„+D • v— r + V D 


lli 


r 2 0R xs 3x 4 x 0 r 2 9x 2 30 2 


uL. u n 9 4 ,i n 3^_ 
r 2 9 x3x 2 90 2 r 4 U 0R 3e 4 


. _2_ n O 

r 2 x0 12 3x 2 3e 2 


(D. 25) 


The equations of motion that result for the shell when the 
matrix (D.25) is used are almost identical to those of Mikulas 
[15], p. 8, or Leissa [16], p. 191* The differences are due to 
assumptions made in the analyses, but the important main diagon- 
al terms are very similar in form. 


To incorporate the effects of the stiffeners in the Peterson 
analysis, the properties of the stiffeners are "smeared-out" , 
l.e., averaged over the shell surface as done in Ref. [15]. An 
equivalent skin thickness, t, is defined by 



where t is the actual skin thickness, A and A„ are the 
s s R 

stringer and ring frame cross-sectional areas, and d and l are 
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the stringer and ring spacings. Using this thickness, the skin 
stiffnesses are computed from 


k 

x 


k e 



5 


and the skin rigidities with 


D„ = D, 


Et 3 

12 ( 1-v 2 ) * 


In (D.25), k xs and k 0R revert back to k x and k 0 (as originally 
used by Peterson) but are computed with the equivalent skin 
thickness. D and D 0R revert back to D x and D 0 everywhere 
except in the third diagonal term and are also computed using 
the equivalent thickness as above. In the third diagonal term D 

X O 

and D 0r are computed from the results in (0.24) where D x and 
D 0 are computed for the actual skin thickness. 


When these modifications are incorporated, it can be seen that 
reasonable approximations are found for the first two diagonal 
terms and for the off-diagonal terms but the skin stiffnesses 
will be somewhat higher than they would have been had k and 
k 0R been used (i.e., the results in (D.24)). 


The Mode Shapes (Eigenvectors) 


The eigenvectors for the fuselage model are the displacement's u, 

v, and w for each eigenvalue, o> . The w component is the 

r — ” 

desired mode shape \p (x) . This mode shape encompasses both the 
plate and shell normal displacements. The maximum normal dis- 
placement computed on either plate or shell is assigned the 
value of 1.0, and all other values of u, v, and w at other 
positions are divided by that maximum value for normalization 
purposes. 
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The mode shapes for the symmetric modes of the shell are 

n* 

^(z,6) = sin C Mn ( "l) n co sn e s (D.26) 

n=u 

and for antisymmetric modes 

n* 

^g(z 3 6) =-sin ^2 C^(-l) n sinn6 . (D.27) 

n=l 

In these two equations, z is again the axial coordinate and 6 is 

Sr 

measured from the bottom of the cylinder. The C Mn are the 
generalized coordinates for the shell (the same as the W Mn of 
Eq . ( D . 7 ) ) • 


The symmetric modes of the floor are 

<(z,x) - sin B« t cJE 

n=0 


cos 


nirx 


P ™* L v Mn ™ L_ 


(D.28) 


and the antisymmetric modes are 


n’ 


*£(«.*> - sin Bp E eft sin Sp 


(D. 29) 


n=0 


Again, z is the axial coordinate measured from the forward end 
of the cylinder and x is the distance measured horizontally in 
the floor plane from the centerline of the fuselage to the 
position of concern on the floor. Lp is the width of the 
floor and is given by 


Lp = 2a sin 0 O > 

pr 

where 0 Q is the floor angle as shown in Fig. 1. The C Mn are the 
generalized coordinates for the floor plate. Note that i|> r (x) is 
the combined set for cylinder and floor, that is. 


^ r (x) = ^(z,e), \p^(z,x) 

& p 


(D. 30) 
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Generalized Mass 


The generalized mass must include all energy in both the 
cylinder and floor and therefore is defined by 


M 


r 



(u r + v r + w r ) m ad6 
s s s s 



(u r + v r + w r ) m dx 
P P P • P 


(D.31) 


where m s and m p are the masses per unit of area of shell and 
plate respectively, a is the cylinder radius, and L is the 
cylinder length. Note that there are two components of the 
modal mass, one for the shell and one for the plate. If a 
particular mode r is, say, predominately a floor mode, will 
be dominated by the second term and the shell contribution will 
then be small. For such a mode, the floor can be seen to re- 
strain the motion induced by fuselage sidewall excitation. 


Sample Output 


Modes are ranked according to the occurrence of their resonance 

frequencies and the values of the generalized coordinates 
sr pr* 

C Mn and CjJ n for each mode are output. For each value of r, 
there is a single value of M that defines the axial mode shape, 
and a sequence of n's that defines the circumferential mode 
shape. Table D-l gives the results for the first 32 modes of the 
1.83m (72 in.) long cylinder used in the present validation study. 
The values of n given in the table are those that contribute most 
to the determination of the mode shape. Ohly 5 terms are retained 
for the shell and 3 for the floor. Fig.D-2 shows an example of a 
typical mode shape, in particular, for the first mode of Table D-l. 
It is apparent from the figure and also from the computed gener- 
ized mass, that the mode shown is basically a floor mode. 
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TABLE D-l. EXAMPLE PROGRAM OUTPUT 


Phase III Test Article: 1.80^n (71 in.) long cylinder, 
stiffened 0.0008m (0.032in) skin, 0.50&n (20 in.) radius with floor at 56.6° 


NODE: 

FREQ 

node 


SHELL 


PLATE 

GENERAL I ZED NASS 

(KG) 

NO 

(HZ) 

TYPE 

h 

N 

CNN 

H 

N 

CNN 

TOTAL 

SHELL W 

PLATE M 

1 

188.52 

SYNN 

1 

2 

-. 06640 

1 

0 

.41496 

1.13889 

.02830 

1. 10253 




1 

3 

-.03362 

1 

2 

.33940 







1 

4 

.02042 

1 

1 

.23432 







1 

1 

.01776 










1 

5 

-.00589 







2 

208.24 

synn 

1 

2 

-.07694 

1 

1 

1.42371 

1.12644 

.05670 

1.05950 




1 

4 

.06646 

1 

2 

-.34033 







1 

5 

-.03793 

1 

0 

-.30912 







1 

3 

-.02065 










1 

1 

-.01360 







3 

217.88 

SYNN 

2 

4 

.01270 

2 

0 

.44544 

1.06241 

.00143 

1.06071 




2 

2 

-.01001 

2 

2 

.41198 







2 

5 

-.00550 

2 

1 

.14239 







2 

1 

.00391 










2 

3 

-.00293 







4 

231.01 

SYNN 

2 

4 

.04328 

2 

1 

1.08434 

1.05432 

.01444 

1.03897 




2 

5 

-.03088 

2 

3 

.20954 







2 

2 

-.01367 

2 

2 

-.13739 







2 

8 

-.00803 










2 

7 

.00740 







5 

293.25 

SYNN 

3 

4 

.00824 

3 

0 

.46520 

1.05532 

.00050 

1.05476 




3 

5 

-.00455 

3 

2 

.43756 







3 

2 

-.00386 

3 

1 

. 10293 







3 

l 

.00195 










3 

7 

.00119 







6 

301.93 

SYNN 

1 

2 

.60434 

1 

2 

.34208 

3.35573 

2.32876 

.47468 




1 

3 

.39210 

1 

1 

.24271 







1 

5 

-.02343 

1 

3 

.11650 







1 

4 

.01026 










1 

6 

.00956 







7 

303.67 

SYNN 

3 

4 

.03058 

3 

1 

1.00245 

1.04019 

.00869 

1.03105 




3 

5 

-.02783 

3 

3 

.19801 







3 

8 

-.00764 

3 

2 

-.08872 







3 

2 

-.00758 










3 

7 

.00698 







6 

318.73 

ANTI 

1 

2 

-.69794 

1 

1 

-.71809 

5.33569 

2.87324 

1.56559 




1 

3 

.33122 

1 

2 

-.37744 







1 

1 

-.17597 

1 

3 

-.13196 







1 

4 

-.10995 










1 

5 

.03684 
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STIFFENED . 

032 IN 

C Yu 

INDER, 20 IN 

RADIUS,. WITH 

FLOOR AT 

56.6 DEGREES 

NODE 

FREQ 

NODE 


SHELL 


PLATE 

GENERALIZED NASS 

(KG) 

NO 

(HZ) 

TYPE 

N 

N 

Cun 

n 

N CNN 

TOTAL 

SHELL H 

PLATE k< 

9 

346. 71 

ANTI 

1 

3 

.65759 

1 

2 -.37430 

3.07613 

2.38329 

.30019 




1 

2 

.29630 

1 

1 -.10529 







1 

1 

.08228 

1 

3 -.02409 







1 

4 

.06592 









1 

5 

-.02470 






10 

434.31 

SYMN 

4 

4 

.00568 

4 

0 .47682 

1.05268 

.00025 

1.05241 




4 

5 

-.00383 

4 

2 .45133 







4 

2 

-.00218 

4 

1 .08093 







4 

1 

.00130 









4 

7 

.00106 






11 

441.78 

SYHH 

4 

5 

-.02558 

4 

1 .96440 

1.03311 

.00632 

1.02654 




4 

4 

.02362 

4 

3 .19308 







4 

8 

-.00753 

4 

2 -.06599 







4 

7 

.00695 









4 

2 

-.00602 






12 

470.53 

ANTI 

2 

3 

.42411 

2 

2 -.48203 

2.27504 

.82788 

1. 34321 




2 

2 

-.05284 

2 

1 -.46964 







2 

4 

-.03870 

2 

3 -.24148 







2 

5 

-.01810 









2 

6 

.01552 






13 

502.23 

SYHH 

1 

1 

-.76763 

1 

2 .66675 

6.31376 

2.76903 

.76399 




1 

4 

-.12362 

1 

0 -.17384 







1 

3 

-.11395 

1 

4 .05223 







1 

6 

.01293 









1 

0 

.01078 






14 

528.78 

ANTI 

1 

4 

-.66794 

1 

1 1.12932 

3.71841 

2.70485 

.23627 




1 

1 

.31690 

1 

2-1.02632 







1 

2 

-.22781 

1 

3 .92727 







1 

5 

.05974 









1 

7 

-.02769 






15 

555.05 

SYMH 

l 

4 

.61552 

1 

3 -.26052 

3.02626 

2.55327 

. 18097 




1 

3 

.37347 

1 

2 -.09653 







1 

2 

-.18507 

1 

1 -.08745 







1 

1 

-.09309 









1 

5 

.07360 






lb 

558.43 

ANTI 

2 

3 

.30944 

2 

1 1.00231 

1 . b5 729 

.48102 

1.10579 




2 

2 

.10406 

2 

3 .96170 







2 

1 

.02332 

2 

2 -.46528 







2 

b 

.01193 









2 

7 

-.00791 
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STIFFENED . 

032 IN 

CYLINOER, 20 IN 

RADIUS. MITH 

FLOOR AT 

56.6 DEGREES 

NODE 

FREQ 

MODE 


SHELL 


PLATE 

GENERALI 2ED HASS 

(KG) 

NO 

( H2 t 

TYPE 

ft 

N 

CNN 

ft 

N 

CMN 

TOTAL 

SHELL tf 

PLATE M 

17 

560.65 

ANTI 

1 

1 

-.25327 

1 

3 

. 55091 

2.00966 

.54454 

1.11528 




1 

5 

.19632 

1 

1 

. 53220 







1 

3 

.11259 

1 

2 

.14949 







1 

2 

-.05185 










1 

5 

.02207 







IS 

572.86 

ANTI 

3 

3 

-.03551 

3 

1 

.55890 

1.22949 

.04255 

1.18227 




3 

5 

.03125 

3 

3 

.46354 







3 

2 

.02289 

3 

2 

.21649 







3 

5 

.01994 










3 

6 

-.01569 







19 

591.88 

SYHM 

2 

3 

.42559 

2 

2 

.06831 

1.73996 

1.52659 

.00197 




2 

2 

.23611 

2 

3 

-.05649 







2 

5 

.27539 

2 

0 

-.02208 







2 

5 

.03191 










2 

6 

-.02568 







20 

607.50 

ANTI 

2 

5 

.44272 

2 

2 

.84374 

2.88669 

1.78872 

.88336 




2 

3 

.32362 

2 

4 

.25752 







2 

2 

.29967 

2 

3 

-.05258 







2 

5 

.06404 










2 

6 

-.05875 







21 

626.01 

ANTI 

3 

3 

-.29317 

3 

2 

1.58856 

1.95333 

.66669 

1.23632 




3 

4 

.20494 

3 

3 

-.54749 







3 

6 

-.09539 

3 

4 

.44915 







3 

5 

.08096 










3 

2 

.04894 







22 

632.02 

ANTI 

1 

1 

.22426 

1 

2 

.64039 

1.99704 

.51331 

1.07277 




1 

2 

-.13579 

1 

4 

.23344 







1 

3 

.12769 

1 

3 

. 20029 







1 

5 

.12610 










1 

4 

.09665 







23 

636.15 

SYHM 

5 

4 

.00462 

5 

0 

. 48042 

1.05173 

.00018 

1.05154 




5 

5 

-.00363 

5 

2 

.45601 







5 

2 

-.00161 

5 

1 

.07378 







5 

1 

.00102 










5 

7 

.00100 







25 

651.20 

SYHH 

5 

5 

-.02553 

5 

1 

.95366 

1.03218 

.00591 

1.02608 




5 

4 

.02158 

5 

3 

. 19080 







5 

8 

-.00763 

5 

2 

-.05895 





5 7 .00715 

5 2 -.00591 
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STIFFENED .032 IN CYLINDER* 20 IN RADIUS, KITH FLOOR AT 56.6 DECREES 


NODE 

FREQ 

NODE 


SHELL 


PLATE 

GENERAL IZEO NASS 

(KG) 

NO 

( H2 ) 

TYPE 

n 

N 

CNN 

N 

N 

CNN 

TOTAL 

SHELL M 

PLATE W 

25 

648.70 

SYNH 

2 

4 

.54137 

2 

3 

-.31842 

3.12167 

2.22892 

.58061 




2 

2 

-.45147 

2 

2 

-. 31722 







2 

5 

.01850 

2 

0 

. 11216 







2 

1 

.01384 










2 

6 

-.00425 







26 

663.30 

ANTI 

4 

3 

-.03615 

4 

1 

.60174 

1. 17010 

.00997 

1.15917 




4 

4 

.01914 

4 

3 

. 53040 







4 

5 

.01472 

4 

2 

.13228 







4 

6 

-.01272 










4 

2 

.01226 







27 

716.23 

ANTI 

4 

4 

.11594 

4 

2 

1.12205 

1.33482 

.16827 

1. 15854 




4 

3 

-.11451 

'4 

4 

.34731 







4 

6 

-.07550 

4 

3 

-. 19925 







4 

5 

.05392 










4 

7 

.03384 







28 

717.33 

ANTI 

2 

2 

-.42010 

2 

2 

.55423 

2.32071 

1.11855 

.98835 




2 

3 

.16269 

2 

4 

.29992 







2 

4 

-.15887 

2 

3 

.25016 







2 

5 

.1 1998 










2 

6 

-.07240 







29 

738. 44 

ANTI 

1 

4 

-.59228 

1 

2 

.81512 

3.63204 

1.85287 

.93833 




1 

3 

.14987 

1 

4 

. 50058 







1 

1 

-.12843 

1 

3 

-.13113 







1 

6 

-.09905 










1 

5 

-.07457 







30 

761.91 

SYNH 

3 

4 

-.63950 

3 

3 

.31731 

2.90922 

2.41279 

.31644 




3 

3 

-.35980 

3 

2 

.15446 







3 

5 

-.00670 

3 

0 

-.05975 







3 

2 

.00615 










3 

1 

-.00485 







31 

779.76 

ANTI 

3 

3 

-.60810 

3 

2 

-.25398 

2.59160 

2.12986 

.26505 




3 

4 

-.29025 

3 

4 

-.16348 







3 

5 

-.10698 

3 

3 

-.16054 







3 

2 

-.08652 










3 

6 

.04237 







32 

811.04 

ANTI 

2 

4 

.44856 

2 

2 

-. 33993 

2.10595 

1.69469 

.20392 




2 

5 

.28040 

2 

4 

-. 28072 







2 

2 

-.25218 

2 

3 

.05882 







2 

1 

-.14779 










2 

3 

-.09074 
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MODE NUMBER 1 
FREQUENCY = 188.51 
SYMMETRIC MODE 
M = 1 
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FIGURE D-2. EXAMPLE SHELL MODE ( Z DEPENDENCY SUPPRESSED) 
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Fig.D-3 shows another mode, the sixth mode of Table D-l. It. is 
found to be basically a shell mode. 



MODE NUMBER 6 
FREQUENCY = 301.93 
SYMMETRIC MODE 
M = 1 



* 

* 

* 


* 

* 

* 

* 

1.0 


* 

* 

* 

* 0.5 

FLOOR * 

DISPLACEMENT 

♦ 

C/L ••••••* • • • .SIDE 


* 


Top C/L 


0.5 

j * 


* 


* 


* 

* 


1.0 


SHELL 

DISPLACEMENT 
(Cylinder Wall) 


* 


* 


* 

. * 

F = * FLOOR LOCATION 
. * 


* 

* 

* 

* 


* 


* 


* 

* 

* 

* 


* * * Bottom C/L 


FIGURE D-3. A SECOND EXAMPLE MODE 


D-19 



APPENDIX E 


MODEL VALIDATION STUDIES 



APPENDIX E - VALIDATION STUDIES 


The models of Section 3 are based on fairly precise solutions of 
the sound transmission problems. Nevertheless, the models will 
have limited prediction accuracy due to errors inherent in the 
input data that is used to describe the various physical systems 
involved, i.e., exterior field, fuselage structure, cabin acous- 
tic space, sidewall trim, etc. The primary goal of the present 
modeling effort has been to minimize these errors through 
elaborate, though limited, system models. Although the overall 
quality that has been achieved is not predictable, it is 
measurable, and a test has been devised for that purpose. 

Test Hardware 


The test configuration is very similar to that represented in 
Figure 1. The actual test rig and hardware are shown in Figures 
E-l through E-3. As seen in Figure E-l, the fuselage and the 
propeller are located downstream of a nozzle that supplies air 
to simulate airplane forward velocity. The fuselage model 
itself is a cylinder 1.83m (72 in.) long and 1.02m (*J0 in.) in 
diameter. The skin is 0.00081m (0.032 in.) thick and is 
stiffened by eighteen (18) stringers spaced on 20° centers. The 
stringers are 90° angles having dimensions of approximately 
0.00953 x 0.0112 x 0.00051m (3/8 x 7/1.6 x 0.020 in.). They are 
riveted to the inside of the skin and pass through cut-outs in 
eight (8) internal ring frames that are spaced along the 
cylinder every 0.2m (8 in.) The frames are aluminum channels 
with dimensions of approximately 0.017 x 0.038 x 0. 00081m (5/8 x 
1-1/2 x 0.032 in. ). 

The cylinder has a structurally integral floor consisting of a 
0. 00081m (0.032 in.) plate stiffened by floor supports of the 
same thickness spaced every 0.2m (8 in.). The supports extend 
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FIGURE E-1. MODEL TEST FACILITY 
(Dimensions in meters) 
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0. 046 m 
(1.8 in.) 


0.00081 m (0.032 in.) skin 
1.016 m (40 in.) diameter 
x 1.83m ( 72 in. ) long 


0.00051 m (0.02 in.) 
stringer, 0.0095 x 
0.011 m (0.375 x 0.44 in.) 


0.039 



0.00081 m (0.032 in.) 
frame, 0.017 x 0.038 n 
(0.67 x 1.5 in.) channel, 
0.2 m (8 in.) spacing. 

8 req'd. 


0,00081 m (0.032 in.) floor beam, 

0.017 x 0.0.038 m (0.67 x 1.5) channel 
1.79 m (70.75 in.) long 


0.00081 m (0.032 in.) floor support at each frame 

0.00081 m (0.032 in.) floor, 1.83 m (72.0 in.) long, 
rivet to supports, bolt to skin, cutout to clear frames 


FIGURE E-2. FUSELAGE MODEL 
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Stiffener 




Clip for end ring and support attachment 


floor 





FIGURE E-3. FLOOR ASSEMBLY 
(Dimensions in meters) 



downward from the floor to the bottom of the cylinder. There 
are also two floor beams (channels of the same dimensions as the 
cylinder ring frames) that run longitudinally, each located 
approximately 0.13m (5.1 in.) from the center of the floor. The 
width of the floor is 0.85m (33*4 in.) leading to a floor angle 
0 q of 56.6 degrees (see Fig.D.l). The outer edge of the floor 
is bolted to the cylinder wall. The cylinder is closed by 
0.013m (1/2 in.) thick end caps that are used to support the 
cylinder in the NASA Langley propeller test facility. The 
entire fuselage assembly is constructed of 2024-T3 aluminum. 

Propeller 


The propeller is a three-b laded, 0.3 scale Hartzell for a Twin 
Otter aircraft with a diameter of 0.76m (30 in.). It is driven 
by a 30 kw (40 horsepower) variable speed electric motor capable 
of turning it up to 8000 rpm. The propeller blades are 
Series 16 airfoils. The geometry of the blades is specified in 
the input data to the ANOPP Propeller Noise Prediction program. 

The data used to define the particular propeller used in the 
present test are proprietary and thus are not included here. In 
the present circumstance the angle-of-attack and local chord 
are specified as a function of radial location. Airfoil coordi- 
nates are specified for several locations and then interpolated 
as required with a cubic spline to fix these variables at all 
locations on the blade. 

Figure E-4 shows the grid coordinates used for the calculations 
required in the ^present test. There are 160 points on the upper 
quarter of the cylinder. In the present test rig, each position 
on the grid, i.e., (k,£) lying in the fuselage surface, has 
coordinates defined by the equivalence relation 

(k,£) ■<->- ( x q , X 2 , x^ ) , 

where (in meters): 
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FIGURE E-4. GRID USED FOR PROPELLER NOISE PREDICTIONS 


- 0.457 + O. 508 . {1-cos [( £-1 )• ir/l8 ] } 
= -0.508 sinC(JUl) •n/18] 

and 

x 1 ^ = 0.622 - 0.089 (k-1) . 


This grid covers all of the upper quarter surface of the 
cylinder forward of the propeller and a somewhat greater amount 
behind it. Because of the lengthy calculations involved in the 
ANOPP program, the data for the lower quarter of the cylinder 
seen by the propeller are obtained from the data for the top 
quarter with the relation (imagining an identical grid below the 
centerline) 

, it it k . n / it it k , x 

p(x 1 ,x 2 , x ^ j t) = p(x 1 ,-x 2 ,x 3> t-T k£ ) , 
where is a time delay given in milliseconds by the result 


kit 


- 333.33 
N 


a 


kit 


N is the propeller rpm and is in degrees and is given by 

the result 


a 


kSL 



The propeller harmonic amplitudes at corresponding points above 
and below the centerline are given by 



bottom 



5 


and the corresponding phases (in degrees) are related by 


kit 

H 


bottom 


= 4 > 


kit 

H 


kit 

+ ~ x H x 360° 
top 1 
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where = BPF is in milliseconds and H is the harmonic index. 
This can also be written as 


<J> 


kH 

H 


= <P 


bottom 


kfc 
H I 


+ 2BHa 


top 


k£ 


where B is the number of propeller blades. 


Conversion to the coordinate system used in Figures 1 and 2 in 
the body of this report is with the relations 


J k 


= Z P ‘ X 3 


and (using 6 = tt/ 2 in Fig. 2) 


\ = I + tan 1 


-x. 


r p " X 1 


The coordinates of the grid point 
equivalence statement 


= <J) p + (£-1) . 

(k,£) are given by the 


(k,Z) (a,e £J z k ) 

In the present test, r p = 0.965m (36 in.) and z p = 0.662m 
(<0*5 in.). The resulting grid has spacing A of approximately 
0.089m (3.5 in.). This spacing is sufficiently close to assure 
a relatively smooth change in phase for each propeller harmonic 
from grid point-to-point. 


Test Description 

The test program discussed in this appendix was conducted at 
NASA Langley Research Center by NASA personnel. The main test 
involved the transmission of propeller noise into the interior 
of the model cylinder. For this test the cylinder and propeller 
were mounted in an anechoic chamber in the configuration shown 
in Figure E-l. Additional tests were performed in a reverbera- 
tion chamber to measure the noise reduction associated with a 
reverberant sound field. Also, decay measurements were made on 
the acoustic field in the cylinder and the structural vibration 
of the cylinder to determine empirical loss factors. 
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The interior of the cylinder was treated with fiber glass blank- 
ets with a density of 9.61 kg/m3 (0.6 lb/ft3). The blankets 
were applied in layers which were nominally 1.3 cm (0.5 inch) 
thick, four layers with a total nominal thickness of 5.1 cm (2 
inch) being applied to the curved surfaces and one layer to the 
end plates. One face of each layer had a vinyl facing with a 
thickness of 0.005 (0.002 inch); the total surface density of 
one layer plus facing was 0.22 kg/m2 (0.045 lb/ft2). On the 
curved surface, three of the fiberglass layers were placed 
between the frames and the fourth layer covered the frame caps. 
The inner surface of the fiberglass treatment on the curved- 
walls was covered with a trim septum consisting of a sheet of 
epoxy/fiberglass NEMA G-10 (0.079 cm or 0.031 inch thick) from 
floor to floor, and a sheet of vinyl (of the same thickness) 
over the upper 120° of the cylinder. The installation is shown 
diagrammatically in Figure E-5. The trim was hard-mounted to 
the floor and attached to the frames by nine soft-mounted 
screws. The total weight of the trim septum was 6.56 kg (14.51 
lb). 

During the propeller noise test, sound levels inside the cylin- 
der were measured using an array of eleven microphones which 
could be located at any selected station along the cylinder and 
could be rotated about the cylinder axis. Measurements were 
made at a total of 196 locations, consisting of 49 locations at 
each of four axial stations. The stations were chosen so that 
the interior of the cylinder was divided into four segments of 
equal volume. Figure E-6 shows the microphone locations on the 
rotating array. The array was positioned at angular locations 
= 0°, + 51.5° and +103°. The radial positions of the micro- 
phones were chosen such that all microphones were associated 
with approximately equal cross-sectional areas. The sound 
levels measured at the different microphone locations were 
averaged on an energy basis to obtain space-average sound 
pressure levels for each measurement axial station and for the 
cylinder as a whole. 
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FIGURE E-5. CROSS-SECTION OF TEST CYLINDER SHOWING TRIM 
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FIGURE E-6. MICROPHONE LOCATIONS IN TEST CYLINDER 
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Measurements were made by NASA at several test conditions, but 
results for only one condition are considered in this appendix. 
The airflow speed into the propeller was 23-d m/s (78 ft/sec) 
and the propeller rotational speed was 4000 rpm. The resulting 
fundamental, or first harmonic, of the blade passage frequency 
is 200 Hz. 

In the case of the noise reduction tests for reverberant excita- 
tion, sound levels were measured at only two axial stations, one 
at the mid-point of the cylinder and the other at the one-sixth 
point. Symmetry was then assumed about the mid-point of the 
cylinder length in order to calculate space-average sound 
pressure levels. 

The bulk of the data reduction was performed in terms of one- 
third octave band spectra. However, a small amount of narrow- 
band analysis was performed in order to obtain a better under- 
standing of the data. 

Measured Interior Sound Levels 


Typical narrowband sound pressure level spectra measured at two 
locations in the cylinder are shown in Figure E- 7 • The two 
selected locations are close to the plane of rotation of the 
propeller (X 3 /D = 0.015 where D is the propeller diameter), 
microphone #1 being near to the center of the cylinder and #8 
near to the periphery (as shown in Figure E-6). Figure E-7 
shows that the contributions associated with the three lowest- 
order harmonics are easily identified but that contributions 
from higher-order harmonics may be difficult to identify. This 
masking of the discrete frequency components by the broadband 
signal limits the upper frequency bound on the useful propeller 
noise data, particularly when, as in the present case, much of 
the data presentation is in terms of one-third octave band 
spectra. 
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Average one-third octave band sound pressure level spectra are 
shown in Figure E-8 for each of the four axial stations of the 
microphone array. Little variation in sound pressure level is 
observed from station to station. The space-average spectrum 
for the cylinder as a whole is plotted in Figure E-9 , which 
also contains the range of measured sound pressure levels. The 
space-average values are compared in Figure E-10 with associated 
95% confidence limits. Because of the large number of samples, 
the 95% confidence intervals are small. 

The spectra show distinct peaks in the frequency bands centered 
at 200, 400 and 630 Hz, the sound levels in these bands being 
controlled by the contributions from the first, second and third 
harmonics of the blade passage frequency. At higher frequen- 
cies, the broadband contributions become important, as indicated 
in Figure E-7, and the one-third octave band spectra are rela- 
tively smooth. The harmonic components may be significantly 
lower in level than are the corresponding one-third octave band 
levels. 

Propeller Noise Field 

The analytical model to predict sound levels in the cylinder 
uses, as data input, a description of the propeller noise field 
(in terms of pressure amplitude and phase) computed for free- 
field conditions using the NASA ANOPP computer program. The 
blocked pressures are then calculated within the present analy- 
tical model using the relationship given in Eq.(43)« 

During the test program at NASA, sound pressure levels were 
measured at several free-field locations surrounding the model 
propeller, and at other locations on a rigid-wall cylinder 
placed close to the propeller. A brief comparison of the 
measured and predicted sound levels has been performed to get an 
indication of the accuracy with which the predictions fit the 
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ONE-THIRD OCTAVE BAND SOUND PRESSURE LEVEL, DB RE 20 MICRO PA 



ONE-THIRD OCTAVE BAND CENTER FREQUENCY, HZ 


FIGURE E-8. AVERAGE SOUND PRESSURE SPECTRA IN CYLINDER 
,AT DIFFERENT MEASUREMENT STATIONS 
(PROPELLER NOISE EXCITATION) 
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ONE-THIRD OCTAVE BAND SOUND PRESSURE LEVEL, DB RE 20 MICRO PR 


100 


90 


80 


70 


60 


50 


40 


30 



> 

'l 
'l 
' I 
' I 
' 1 
1 1 


SPACE-A 
MAX. AN 

i'ERACE S 
3 MIN. f- 

PL 

EASURED 

VALUES 


1 1 
' A > 

‘ \ > 
'i\' 

• / \ i 

'Mi 

' / \ > 
1 l » 
' / 1 > 

'i 
'i 
' i 
' i 
' i 
' i 








/’ 

/ ; 

' \ 





■ 

■ 

■ 





’ \ 
l \ 

/ \ 

* \ 

/ t 

/ l 

/ \ 


■ 

1 










1 1 

1 1 

1 1 

, , 

1 

\ 

\ 

\ 

\ 

1 l 

i i 

i i 


63 125 250 500 1000 2000 4000 8000 

ONE-THIRD OCTAVE BAND CENTER FREQUENCY, HZ 


FIGURE E-9. SPACE-AVERAGE AND RANGE OF VALUES FOR 
PROPELLER-INDUCED SOUND LEVELS INSIDE 
TEST CYLINDER 
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ONE-THIRD OCTAVE BAND SOUND PRESSURE LEVEL, DB RE 20 MICRO PA 
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FIGURE E-10. SPACE-AVERAGE AND 95% CONFIDENCE LIMITS FOR 
PROPELLER-INDUCED SOUND LEVELS INSIDE 
TEST CYLINDER 
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test environment. This comparison was made for axial locations 
associated with grid line £ = 1 shown in Figure E-4. The com- 
parison was restricted to pressure level, because no phase data 
were available from'the test program. 

The comparison of measured and predicted free-field pressure 
levels is shown in Figure E-ll. In this case both the measure- 
ments and the predictions were performed by NASA personnel. The 
agreement is generally good except at large distances from the 
plane of rotation where the pressure levels are relatively low 
and could be affected by broadband flow noise. 

Figure E-12 contains the corresponding comparison for the 
pressure field on a rigid cylinder (i.e.,the blocked pressure). 
The test data were obtained directly from NASA measurements and 
the predictions from the application of Eq.(43) to the output of 
ANOPP. The equation provides an empirical relationship for 
calculating the effect of pressure reflections at tne surface of 
the cylinder. In general the predicted levels are higher than 
the measured values, the differences being 0 to 5 dB at loca- 
tions in the neighborhood of the plane of rotation where the 
pressure levels are highest. 

Acoustic Loss Factors 


Acoustic reverberation decay measurements were made inside the 
test cylinder using pink noise and sinusoidal excitations. 
Acoustic absorption coefficients were computed from the 
reverberation time Tr using the relationship 


a 


17.6ttV 

c ST_ 
o R 


where S and V are,., respectively, the surface area and volume of 
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FIGURE E-ll. MEASURED AND PREDICTED FREE-FIELD 
SOUND LEVELS FOR TEST PROPELLER 
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FIGURE E-12. MEASURED AND PREDICTED BLOCKED SOUND LEVELS 
INDUCED BY PROPELLER ON TEST CYLINDER 






the interior of the cylinder. Values of a, averaged over 
approximate one-third octave bands, are given in Figure E-13- 

The acoustic loss factor n n was computed from the equation 

c Sa 
n = o 

'n Ta y\T 

Resulting values for the acoustic loss factor are plotted in 
Figure E-l^. 

Also contained in Figure E-14 are band-average values of the 
acoustic loss factor calculated using the analytical model out- 
lined in Appendix A. It is necessary in this analytical model 
to assign a value to the structural loss factor nip of the trim 
septum. A' value of 0.5 was assumed in [7] for the trim loss 
factor, on the basis of the best fit with test data then avail- 
able. The same value was assumed initially for the present 
analysis. At frequencies above 500 Hz, the measured data are 
scattered about the predicted spectrum, but at lower frequencies 
the predicted values are significantly higher than the measure- 
ments . 

When the trim structural loss factor was increased by a factor 
of 2 to a value of 1.0, the calculated acoustic loss factor 
decreased in the 250 Hz band (which contains a predicted reson- 
ance frequency of the trim) but the change at other frequencies 
was small or negligible. The disagreement between measurements 
and predictions in the lower frequency -range was not affected to 
any significant extent. 

Structural Loss Factor 


The damping of the cylinder structure was determined experi- 
mentally from vibration decay measurements using sinusoidal 
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FIGURE E-1 3. MEASURED ACOUSTIC ABSORPTION COEFFICIENTS 
IN TEST CYLINDER 
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excitation. The resulting empirical values are plotted in 
Figure E.15, where they are compared with band-average values 
predicted on the basis of the analytical model given in 
Appendix A. The analytical model includes the influence of the 
sidewall trim and requires, as an input, an assigned value for 
the trim structural loss factor n T . The two calculated loss 
factor spectra in Figure E-15 assumed trim loss factors of 0.5 
and 1.0. 

The comparison of measured and predicted values shows good 
agreement at frequencies above 400 Hz, but the single experimen- 
tal data point at lower frequencies is an order of magnitude 
smaller than the corresponding predicted values. For diagnostic 
purposes, an alternative analytical representation was construc- 
ted whereby the structural loss factor was arbitrarily limited 
to a maximum value of 0.15 in the frequency band below 400 Hz. 

Predicted Interior Sound Levels 


The model scale experiments discussed in the preceding sections 
of this appendix were performed in order to obtain data for 
comparison with sound levels predicted using the analytical 
model described in this report. The predictions were made 
using a computer program based on the analytical model outlined 
in this report. 

Initially, the interior sound levels were calculated without 
placing any restraints on the structural loss factor for the 
cylinder and with an assumed structural loss factor for the trim 
of n T = 0.5. Space-average sound pressure levels were computed 
for the five lowest-order harmonics of the propeller blade 
passage frequency. These levels are compared in Figures E-16 
with the range of measured values and the associated space- 
average levels, and in Figure E-17 with the average and 95% 
confidence limits. The predicted sound levels refer to discrete 
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FIGURE E-1 6. COMPARISON OF PREDICTED AND MEASURED SOUND 
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FIGURE E-17. COMPARISON OF PREDICTED AND MEASURED SOUND 
1EVELS IN CYLINDER INDUCED BY PROPELLER 
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frequencies whereas the measured levels are associated with 
one-third octave bands. However the harmonic frequencies of 
200, 400, 600, 800 and 1000 Hz are centered within the corre- 
sponding one-third octave bands so that there should be no 
filter cut-off problems. 

The comparisons show that the predicted levels are generally 
lower than the corresponding measured one-third octave band 
space-average levels. The largest difference occurs in the 
highest order harmonic considered (1000 Hz) where the measured 
levels may be dominated by broadband noise rather than the dis- 
crete frequency component associated with the propeller blade 
passage sound. 

A small number of parametric studies have been performed in 
order to assess the sensitivity of the analytical model to 
various factors. These studies were limited to considerations 
of the changes in trim and cylinder structural loss factors 
discussed previously. In summary, the trim structural loss 
factor was assigned a value of 1.0 and the cylinder structural 
loss factor was limited to a maximum value of 0.15. The result- 
ing predicted interior sound levels are shown in Figure E-18. 

The effect of the changes of the loss factors was small except 
for the first harmonic at 200 Hz. The predicted level at this 
frequency now lies within 2 dB of the measured one-third octave 
band space-average value, instead of being 10 dB lower as shown 
in Figure E-17 for the initial calculations. 

Further diagnostic analysis is highly desirable, particularly 
with respect to the model for the sidewall treatment, but such 
analysis was not possible in the present study. 
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FIGURE E- 18. COMPARISON OF PREDICTED AND MEASURED SOUND 
LEVELS IN CYLINDER INDUCED BY PROPELLER (TRIM 
LOSS FACTOR = 1.0, TOTAL STRUCTURAL LOSS 
FACTOR LIMITED TO 0.15 MAXIMUM) 


E-29 







Noise Reduction 


In addition to predicting the transmission of propeller noise 
into an airplane fuselage, the analytical model also contains a 
capability for calculating the noise reduction associated with 
a reverberant incident sound field. This capability has been 
the subject of validation studies in earlier phases of the 
development of the analytical model [3,7]. Additional compari- 
sons are discussed here. 

Interior sound pressure levels were measured at two axial 
stations inside the cylinder, and average noise reductions 
computed for the two stations. The resulting noise reduction 
spectra are plotted in Figures E-19 and E-20* At frequencies of 
160 Hz and above the two average spectra have similar values, as 
was the case in Figure E-8 for propeller noise excitation. 

There are, however, large differences in noise reductions for 
the two stations at frequencies of 80 and 100 Hz. 

Space-average noise reductions for the cylinder as a whole were 
computed assuming that the sound field in the cylinder was 
symmetrical about the mid-point of the cylinder length. The 
space-average spectrum is shown in Figure E-21 with the range of 
measured values and in Figure E-22 with 95% confidence limits. 

Figures E-21 and E-22 also contain the space-average noise 
reductions which were computed. It was assumed for these calcu- 
lations that the trim structural loss factor was n T = 1.0 and 
the structural loss factor was limited to a maximum value of 
0.15. In addition, allowance was made for noise transmission 
through the end plates; the method used was the same as that in 
[7] and it was applied to the four lowest frequency bands, 80 to 
160 Hz. Without this modification the analytical model would 
predict significantly higher noise reductions at low frequen- 
cies . 
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FJGURE E-19. AVERAGE NOISE REDUCTION MEASURED AT STATION 
AT 50% OF CYLINDER LENGTH 
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FIGURE E-20. AVERAGE NOISE REDUCTION MEASURED AT STATION 
AT 83% OF CYLINDER LENGTH 
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FIGURE E-21. COMPARISON OF MEASURED AND PREDICTED SPACE- 
AVERAGE NOISE REDUCTION FOR TEST CYLINDER 
(AVERAGE AND RANGE OF VALUES , TRIM LOSS 
FACTOR = 1.0. TOTAL STRUCTURAL LOSS FACTOR 
LIMITED TO 0.15 MAXIMUM) 
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FIGURE E-22. COMPARISON OF MEASURED AND PREDICTED SPACE- 
AVERAGE NOISE REDUCTION FOR TEST CYLINDER 
(AVERAGE AND 95% CONFIDENCE LIMITS, TRIM LOSS 
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The computer program assumes that no noise is transmitted 
through the ends of the cylinder, since such contributions are 
not usually important in aircraft. Furthermore, because of the 
directivity characteristics of propeller noise, the noise levels 
on the end plates of the test cylinder would be low when the 
cylinder was exposed to the propeller excitation. Thus, noise 
transmission through the end plates would be negligible in the 
propeller noise test. The contribution becomes significant 
only for reverberant field excitation. 

The agreement between measured and predicted noise reductions is 
not as good as was the case for propeller noise excitation. The 
discrepancy between predicted and measured noise reductions at 
high frequencies was observed in previous data [7] where it was 
believed to be due to flanking paths through exposed stiffeners. 
There are no flanking paths of this type in the present test 
model. Two possible explanations for the discrepancy are (a) 
other flanking paths are present, and (b) the analytical model 
over-estimates the influence of the sidewall treatment. 

The noise transmission measurements were repeated with the 
fiberglass treatment on the end plates increased from a thick- 
ness of 1.3 cm (0.5 inch) to 3.1 cm (2 inches). The resulting 
change in measured space-average noise reduction is shown ir 
Figure E-23 which contains spectra associated with the two 
treatments. Figure E-23 contains corresponding predicted noise 
reduction spectra. The two figures show somewhat different 
trends with the predictions giving an increase in noise reduc- 
tion in the frequency range 315 to 1250 Hz whereas the increase 
in measured noise reduction occurs mainly at frequencies above 
800 Hz. 

The predicted spectra in Figure E-24 are influenced directly by 
the analytical model for the fiber glass material on the end 
plates. Figure E-25 compares the acoustic loss factors computed 
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FIGURE E-23. MEASURED NOISE REDUCTIONS FOR CYLINDER WITH 
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for the interior of the cylinder for the two end plate treat- 
ments. The results show that a resonance-type effect predicted 
for the trim is shifted down in frequency from 1600 Hz for the 
1.3 cm treatment on the end plates to 600 Hz for the 5*1 cm 
treatment . 

General Comments 


This brief discussion of validation studies for the analytical 
model has shown fairly good agreement between predicted and 
measured sound levels in a cylinder exposed to propeller noise 
excitation. The discussion has also identified several items 
for future study in order to improve the analytical model. For 
the propeller noise pressure field, the procedure for estimating 
pressure reflection effects should be reviewed. Also, the ana- 
lytical model for relative phase between the upper and lower 
grids should be validated. 

The present analysis represents the first time that the sidewall 
treatment has been incorprated as an integral part of the sound 
transmission model. As such, the analysis makes a major step in 
providing a comprehensive model for airplane interior noise. 

The validation studies indicate that some improvements to the 
sidewall treatment model are required in order to improve the 
agreement between predictions and measurements. It is recom- 
mended that these improvements be made. 
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LIST OF SYMBOLS 



LIST OF SYMBOLS 


A Cylinder surface area, used in Eq.(4). 

or 

A Interior (cylinder & floor) surface area, used in 

conjunction with structural/acoustic coupling 
function f ' (n, r ) 

A' Transmitting area of cylinder without trim 

A Exterior cylinder surface area 

A Transmitting area of cylinder with trim 

t 

(A +A') Total transmitting area of cylinder, floor to floor 

A™ Amplitude of Fourier component of blocked propeller 

pressure signature at propeller harmonic H and grid 
location m = (k,£), see Eq.(43) 

A m Amplitude of Fourier component of free field propel- 

ler pressure at harmonic H and grid location m=(k,£) 

a Radius of cylinder 

a™ Fourier series coeffficient of propeller pressure, 

H 

for harmonic H at location m; see Eqs . (27 ),(29) 

Mean propeller pressure amplitude at location m 
(defined in Eq.(28)) 

Functions defined in Section 3*3, Eq.(17), et seq 
Number of propeller blades 

N33 

BPF Propeller blade passage frequency (Hz); BPF = 
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LIST OF SYMBOLS 

(Continued ) 

Functions defined in Section 3*3, Eq.(l?) et seq 

Fourier series coefficient of propeller pressure 
for harmonic H at location m; see Eqs . (27) , (29) 

Function defined in Eq. (76) 

Cosine integral; Cin(z) =f (l-cost)dt/t 

J o 

Floor and shell generalized coordinates for 
structural mode r = (M,N); see Eqs. (46)-(49) 

Cpbi (^1 * ’ ) Cospectral density function of the blocked 

exterior pressure field 

Cp(x| x';o)) Cospectral density function of the exterior 
pressure field 

R I 

C W ,C W ,C W Trim parameter, derived from the trim transfer 
matrix, Eq.(A.8); C w = C w + iC^ 

C (£,co)C (c,w) Cospectral density functions of the exterior 

x y 

pressure field in the axial and transverse 
directions respectively; see Eq.(60) 

c Q Speed of sound in air 

c n , c r Functions defined in Section 3*3, Eq.(17) et seq 

Constant percentage bandwidth parameter, where 

Aw = c a) [c = 0.232 for one-third octave bands] 
co w 


6 n > b r 
b H 

C( J) 
Cin( z ) 


pr sr 
°Mn°Mn 
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LIST OP SYMBOLS 

( Continued ) 

D njo Function defined in Section 3*3, Eq.(17) et seq 

E[ ] Expected value of a function 

f'(n,r) Interior structural/acoustic coupling factor; see 

Eq. (45) 

f ' (n, r)=f ' (qi, r) Interior structural/acoustic coupling factor 

including effect of trim factor see 
Eq. (52) 

f^ Frequency of propeller 1st harmonic; f^ = 1 /T q = 

BPF 

f H Frequency of propeller harmonic H; f H = Hf.^ 

fq m Acoustic/structural coupling factor in axial 

direction; see Eq.(57) 

g n ,g r Functions defined in Section 3.2, page 3-11 

H Propeller harmonic order, used as superscript to 

denote functions evaluated at frequency w 

II 

i Acoustic mode number counter for fuselage cross- 

section modes, associated with mode n =(q,i) 

I 1 I 2 I^ Integrals defined in Equations (63) and (70)— (73) 

j Circumferential location on fuselage wall, 8., a 

J 

boundary point at which the acoustic eigenvector 
is evaluated (see Fig.C-2) 
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LIST OP SYMBOLS 

(Continued) 



Structural joint acceptance function in axial 
direction 



Structural joint acceptance function in 
circumferential direction 


j* (w) Structural joint acceptance in axial and 

circumferential directions; j£(w) = j Apj ( 03 ) 
= see Eq.(5B) 

2 

A rev 

j (co) Joint acceptance for reverberant/diffuse 
excitation 


<j 2 (w)> 


rev 

r 


k 

or 

k 


L 

L 

P 

H 


In , In 
n 5 r 


M 


Joint acceptance for reverberant excitation 
averaged over structural modes resonant in band Aw 

Acoustic wave number, k = 2 tt/A 

Axial non-dimensional coordinate for grid point; 
see Figure E— ^ 

Fuselage structure length 

Floor width (wall to wall) 

Circumferential non-dimensional coordinate for 
grid point; see Figure E-*)) 

Functions defined in Section 3«3, Eq.17 et seq 

Number of axial half-wavelengths for structural 
mode r = (M,N) 



LIST OP SYMBOLS 

(Continued ) 

Generalized modal mass, for structure mode r 

m = (k,£) Grid point on surface of cylinder used for 
propeller noise predictions; see Figure E-4 
or 

m Average surface mass/unit area of cylinder 

N Structural mode counter, associated with mode 

r = (M,N ) 
or 

N Propeller rpm 

N n ,N r Number of acoustic modes or structural modes in 

frequency band Am 

n Symbolizes acoustic mode n = (q,i) 

or 

n Number of circumferential wavelengths (or 

transverse half-wavelengths) in fuselage 
shell (or floor); see Eqs. 46-49. 

n* Number of terms in displacement series for 

fuselage shell (or floor) 

n b ,n^ Number of boundary points on the fuselage shell 

(or floor) at which the acoustic eigenvectors are 
defined 

n n Modal density of acoustic modes 

n r Modal density of structural modes 
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P(«n) 

P ( % } 

P b l ( x, t) 


<Pbl> 

<PiU,u)> t 


<p? > 

1 S , t 


< ^6 >5 , t 


<p 2 > 

n Sjt 



q 


r 


r 

P 


LIST OP SYMBOLS 

(Continued) 

Probability distribution function for to in Aw 
J n 

Probability density function for in Aw 

Exterior pressure over the blocked (immobile) 
fuselage 

Band-limited mean square blocked pressure 

Interior mean square pressure at location £ 

Space-averaged band-limited mean square interior 
pressure 

Space-averaged band-limited mean square exterior 
pressure for a reverberant field 

Space-averaged band-limited mean square modal 
pressure, for nth mode in interior volume V 

Function defined in Section 3.2, Eq.(12) 

Number of axial half-wavelengths for acoustic mode 
n a (q , i) 

Symbolizes structural mode r a (M,N) 

Radial distance from center of fuselage cylinder 
to the axis of rotation of the propeller. 


F-6 



LIST OP SYMBOLS 

(Continued ) 

W* 1 *' ;u>) Average cross correlation of the exterior 
blocked pressure over the fuselage 

S Absorbing surface area of fuselage sidewall 

S e Absorbing surface area on each end surface 

(bulkhead ) 

Sp(u>) Power spectral density of exterior pressure 

Spki(w) Power spectral density of exterior blocked 

pressure 

S pb i (x | x ' ;oo) Cross spectral density of exterior blocked 
pressure 

Si(z) Sine integral; Si(z) = / z sint.dt/t 

o 

t time 

T Period of rotation of propeller; T = 60/N 

T o ,Tp Period of propeller noise signature; T q = = 

(BPF) -1 = T/B 

u In-plane axial displacement of cylinder wall (or 

floor) 

V Volume of cavity 

v Circumferential (or transverse) displacement of 

cylinder wall (or floor) 

W diss Power dissipated on the cabin walls 
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w. 


in 


W 


int 

rad 


(to) 


W 


int 

abs 


(to) 


w 


X 


X 


m 


x^X2X^ 


LIST OP SYMBOLS 

(Continued ) 

Net power inflow 

Spectral density of power radiated by structure 
into interior acoustic space 

Spectral density of power absorbed on inner wall 
of the space from interior acoustic field 

Cylinder wall (or floor) normal displacement 

Transverse coordinate; see Figure 1 

Location on exterior surface of fuselage 

Location of grid point on exterior surface of 
fuselage 

Local coordinate systems; see Figure E-4 


y Vertical coordinate, relative to fuselage 

centerline (see Figure 1) 

z Axial coordinate, relative to forward end of the 

fuselage structure (of length L); see Figure E-4 

Axial coordinate for grid point k, see Figure E-4 

z Location of propeller relative to the forward end 

of the fuselage structure (of length L); see 
Figure 1) 
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LIST OF SYMBOLS 

(Continued ) 

Band average absorption coefficient 

Incidence angle between propeller and location x^; 
see Figure 2 (in degrees) 

Grid spacing for propeller noise predictions 

Area associated with each grid point; see Figure E-4 


Aoi( radians/sec) Frequency band of width Aw = c^w 

n<Aw symbolizes modes resonant below band 
neAw symbolizes modes resonant inside band 
n>Am symbolizes modes resonant above band 


S( ) 

£ n 


Q 


n 


n 


n 


r 


t 



n 


T ! 

r 


Delta function 

V// v 4)^dv Acoustic mode normalization factor 

Acoustic mode normalizing factor in axial direction 
(see Eq . C . 1 1 ) 

Transverse coordinate; see Section 3*6 
Acoustic mode loss factor 
Structural mode loss factor 

Structural loss factor, including damping due to 
trim; Eq . (82) 

Internal radiation loss factor, due to closely 
coupled structural and acoustic modes; Eq.(83) 
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LIST OP SYMBOLS 
(Continued ) 


n n Average one-third octave band acoustic mode loss 

factor 


h = n'+n' ,+ n®^ Average one-third octave band structural mode 

I ± Jl 1/a.Q 

loss factor 

^struc Average one-third octave band structural loss 

factor 


rad 


Average one-third octave band radiation loss factor 


int ext 
n rad 5 n rad 


Average one-third octave band internal and external 
radiation loss factors defined after Eq. (19) 


e 


Angular coordinate, relative to fuselage bottom 
centerline; see Figure 1 


e 


£ 


Angular coordinate for grid location (k,£ ) 




Angle at which fuselage shell/floor joint is 
located 


6 


ir° 




K. 

or 

K 


Angle 0 for point j on fuselage wall, a boundary 
point for the acoustic eigenvectors 

Angles defining mid-points between boundary point j 
and adjacent boundary points 

Axial coordinate; see Section 3*6 

Conductance for trim on end surface of cylinder 
interior 
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ffi 3 


LIST OF SYMBOLS 

(Continued ) 


| Interior cavity location 

£ e Conductance for trim on cylinder (fuselage) sidewall 

p Density of air inside the cylinder 


Density of air outside the cylinder 


or 


Time delay for cross-correlation 


Acoustic transmission coefficient for diffuse field 

excitation; x = x„ + x 

r R 

Field incidence transmission coefficient for mass 
controlled panels; defined in Eq. (20) 


'mL 


\ mto J 


T t 


Mass law sound transmission coefficient 
Trim transmission coefficient, defined in Eq.(A.22) 


R 


V 1 ’ 


Resonance transmission coefficient for diffuse 
field, defined in Eq.(19) 

Generalized mass for two-dimensional acoustic mode 
i, defined in Appendix C 


<f> Angular position of propeller hub relative to 

fuselage bottom centerline; see Figure 2 

Phase of Fourier component of propeller pressure 
signature at propeller harmonic H and grid location 
m3 (k, l) 
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LIST OF SYMBOLS 

(Continued) 

<f> n (C) = 4>q ^ ( C ) Mode shape, or eigenfunction, of the nth mode 

cf the cavity at location £ 

4)^(0^) Mode shape of ith acoustic mode of the fuselage 

cross-section evaluated on the fuselage wall at 
location j, angle 6. 

J 

Generalized modal forcing function due to propel- 
ler noise, mode r at propeller harmonic H; see 
Section 3.4 

<4 (x) Mode shape, or eigenfunction, of the rth mode of the 

structure, at location x 

4>p(z,x) Floor displacement in structure mode r 

ip^(z,6) Fuselage shell displacement in structure mode r 

s 

Angular frequency (rads/sec) 

03 Angular frequency of propeller harmonic H 

H 

u) n Acoustic mode resonance angular frequency 

03 ^ Structure mode resonance angular frequency 

< > Band-limited, space-averaged and time-averaged value 
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